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Stratospheric  horizontal  wavenumber  spectra  of  winds, 
potential  temperature,  and  atmospheric  tracers 
observed  by  high-altitude  aircraft 
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Abstract.  Horizontal  wavenumber  power  spectra  of  vertical  and  horizontal  wind  velocities, 
potential  temperatures,  and  ozone  and  N2O  mixing  ratios,  as  measured  in  the  mid-stratosphere 
during  73  ER-2  flights  (altitude  ss  20km)  are  presented.  The  velocity  and  potential  temperature 
spectra  in  the  100  to  1-km  wavelength  range  deviate  significantly  from  the  uniform  —5/3  power 
law  expected  for  the  inverse  energy-cascade  regime  of  two-dimensional  turbulence  and  also  for 
inertial-range,  three-dimensional  turbulence.  Instead,  steeper  spectra  approximately  consistent 
with  a  -3  power  law  are  observed  at  horizontal  scales  smaller  than  3  km  for  all  velocity  components 
as  well  as  potential  temperature.  Shallower  spectra  are  observed  at  scales  longer  than  6  km.  For 
horizontal  velocity  and  potential  temperature  the  spectral  indices  at  longer  scales  are  between 
-1.5  and  -2.0.  For  vertical  velocity  the  spectrum  at  longer  scales  becomes  flat.  It  is  argued  that 
the  observed  velocity  and  potential  temperature  spectra  are  consistent  with  gravity  waves.  At 
smaller  scales,  the  shapes  are  also  superficially  consistent  with  a  Lumley-Shur-Weinstock  buoyant 
subrange  of  turbulence  and/or  nonlinear  gravity  waves.  Contemporaneous  spectra  of  ozone  and 
N2O  mixing  ratio  in  the  100  to  1-km  wavelength  range  do  conform  to  an  approximately  uniform 
—5/3  power  law.  It  is  argued  that  this  may  reflect  interactions  between  gravity  wave  air-parcel 
displacements  and  laminar  or  filamentary  structures  in  the  trace  gas  mixing  ratio  field  produced 
by  enstrophy-cascading  two-dimensional  turbulence. 


1.  Introduction 

Studies  of  data  collected  by  commercial  aircraft  fly¬ 
ing  in  the  troposphere  and  lower  stratosphere  (altitudes 
:$14  km)  during  the  Global  Atmospheric  Sampling  Pro¬ 
gram  (GASP)  revealed  a  —5/3  power  law  for  the  power 
spectrum  of  horizontal  velocity  as  a  function  of  horizon¬ 
tal  wavenumber  over  a  wide  range  of  scales  [Gage,  1979; 
Lilly  and  Petersen,  1983;  Nastrom  et  al.,  1984;  Nastrom 
and  Gage,  1985].  Although  this  spectral  shape  is  consis¬ 
tent  with  the  inertial  range  of  three-dimensional  (3-D) 
homogeneous  turbulence  [Tennekes  and  Lumley,  1972], 
the  scales  at  which  these  spectral  shapes  are  observed 
(up  to  1000  km)  are  far  longer  than  can  conceivably  be 
associated  with  3-D  turbulence  in  the  atmosphere. 

Gage  [1979]  and  Lilly  [1983]  proposed  that  these  ob¬ 
servations  reflect  upscale  transfer  of  turbulent  energy 
that  has  been  “two-dimensionalized”  by  atmospheric 
stratification.  This  interpretation  was  based,  at  least 
in  part,  on  laboratory  studies  of  stratified  turbulence 
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[e.g.,  Lin  and  Pao,  1979],  which  showed  that  initially- 
homogenous,  3-D  turbulence  in  a  stratified  tank  was 
rapidly  transformed  into  quasi-2-D  motion  with  pre¬ 
dominantly  vertical  vorticity.  In  2-D  turbulence,  ki¬ 
netic  energy  injected  at  a  certain  scale  will  cascade  to 
larger  scales  while  enstrophy  cascades  to  smaller  scales 
[Kraichnan,  1967].  This  results  in  a  kinetic  energy  spec¬ 
trum  obeying  a  —5/3  power  law  for  scales  longer  than 
the  scale  of  the  injection,  while  a  —3  power  law  holds 
for  scales  smaller  than  the  injection  scale.  In  the  picture 
proposed  by  Gage  and  Lilly,  energy  is  injected  into  2-D 
mesoscale  turbulence  by  even  smaller-scale  3-D  motions 
(e.g.  convective  bubbles)  which  have  been  flattened  by 
stratification  forces.  A  similar  process  is  thought  to 
operate  at  much  larger  scales  where  energy  is  injected 
into  2-D  turbulence  at  scales  ~5000  km  by  baroclinic 
waves  [Charney,  1971].  This  results  in  a  kinetic  en¬ 
ergy  spectrum  obeying  a  —3  power  law  at  scales  shorter 
than  ~5000  km.  Given  continuous  energy  input,  the 
upward-cascading  energy  from  small-scale  energy  injec¬ 
tions  proposed  by  Gage  and  Lilly  should  encounter  the 
downward  cascading  enstrophy  from  the  large-scale  in¬ 
puts  at  some  intermediate  length  scale.  Spectra  deter¬ 
mined  from  the  longest  flights  in  the  GASP  data  do  in 
fact  show  a  change  in  slope  from  —3  to  —5/3  at  scales 
~1000  km  [Nastrom  and  Gage,  1985;  Gage  and  Nas¬ 
trom,  1985].  Horizontal  wavenumber  spectra  for  verti¬ 
cal  velocities  were  not  presented  by  the  GASP  investi¬ 
gators.  Theoretical  predictions  for  the  vertical  velocity 
spectrum  based  on  2-D  mesoscale  turbulence  theories 
are  not  available. 

Earlier  studies  found  steeper  spectral  indices  for  hor¬ 
izontal  velocities  at  smaller  horizontal  scales  within 
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strongly  stratified  environments,  such  as  above  a  low- 
level  inversion  [Koprov,  1965;  Myrup,  1968;  Berman, 
1976]  and  in  the  upper  troposphere  and  lower  strato¬ 
sphere  [Shur,  1962,  1971;  Vinnichenko  et  al,  1965; 
Lilly  and  Lester,  1974],  These  studies  reported  spec¬ 
tra  with  indices  closer  to  —3  than  —5/3  when  strong 
stratification  was  present.  Nastrom  et  al.  [1987]  found 
intensification  and  similar  steepening  of  GASP  spec¬ 
tra  at  large  wavenumbers  over  mountainous  terrain. 
Horizontal  wavenumber  spectra  of  vertical  velocities 
were  also  presented  in  some  of  these  studies  [ Koprov , 
1965;  Myrup,  1968;  Shur,  1971]  and  were  found  to  have 
spectral  indices  at  small  scales  similar  to  those  of  the 
horizontal-velocity  spectra.  The  explanation  proposed 
for  these  findings  was  that  strong  stratification  led  to 
significant  loss  of  energy  as  turbulent  motions  worked 
against  buoyancy  forces.  This  “buoyant  subrange”  of 
turbulence  was  proposed  initially  by  Bolgiano  [1962] 
,Shur  [1962],  and  Lumley  [1964],  Power  law  spectra 
~jfcr  were  predicted,  where  k  is  the  scalar  wavenum¬ 
ber,  and  anticipated  spectral  indices  r  ranged  between 
—  11/5  [Bolgiano,  1962]  and  —3  [Lumley,  1964],  with 
the  latter  value  more  widely  cited.  Weinstock  [1978] 
later  argued  that  a  universal  slope  does  not  exist  for 
buoyancy-modified  turbulence,  and  derived  a  complex 
expression  for  the  spectral  energy  density  which  de¬ 
pends  on  the  flux  Richardson  number  and  the  local  rms 
velocity  fluctuations.  A  wider  range  of  spectral  indices, 
— 4/3  >  r  >  — 3,  was  obtained  than  in  the  earlier  stud¬ 
ies.  Weinstock  [1985]  and  Dalaudier  and  Sidi  [1987]  de¬ 
veloped  the  theory  of  potential  energy  and  temperature 
spectra  in  the  buoyant  subrange  and  predicted  spectral 
indices  around  -3. 

At  intermediate  to  small  scales  (1-100  km),  mo¬ 
tions  which  are  thought  to  be  gravity  waves  are  fre¬ 
quently  detected  in  stratospheric  aircraft  data  [e.g., 
Lilly  and  Lester,  1974;  Nastrom  et  al.,  1987;  Gary, 


1989;  Bacmeister  et  al.,  1990ab,  1994;  Danielsen  et  al., 
1991;  Nastrom  and  Fritts,  1992;  Pfister  et  al.,  1993; 
Chan  et  al.,  1993;  Alexander  and  Pfister,  1995].  Sev¬ 
eral  of  these  studies  examined  statistical  enhancements 
in  variance  over  orography  detected  by  aircraft  [Nas¬ 
trom  et  al.,  1987;  Bacmeister  et  al.,  1990b;  Nastrom 
and  Fritts,  1992].  The  remaining  studies  examined  in¬ 
dividual  wave  events  associated  with  clearly  identifiable 
forcing  from  orography  or  convection.  Spectral  analy¬ 
ses  of  airborne  lidar  measurements  of  mesospheric  vari¬ 
ability  have  been  found  to  be  consistent  with  a  gravity 
wave  interpretation  [e.g.,  Fritts  et  al.,  1989;  Kwon  et  al, 
1990;  Hostetler  and  Gardner,  1994],  However,  the  rela¬ 
tive  contributions  of  gravity  waves,  large-scale  motions, 
and  2-D  and  3-D  turbulence  to  stratospheric  horizontal 
wavenumber  spectra  are  still  not  clear  [e.g.,  Nastrom 
et  al. ,  1987] .  In  this  study  we  will  attempt  to  to  charac¬ 
terize  the  climatological  structure  of  variablility  in  the 
lower  stratosphere  using  high-resolution  meteorological 
and  trace  gas  measurements.  We  will  attempt  to  deter¬ 
mine  the  importance  of  gravity  wave  motions  versus  2-D 
and  3-D  turbulent  motions  in  producing  this  variability. 

2.  Data  and  Methods 

2.1  Description  of  Data 

In  this  study  we  will  compute  spectra  of  horizon¬ 
tal  and  vertical  velocities,  potential  temperatures,  and 
ozone  and  N2O  mixing  ratios,  using  data  obtained  in  the 
midstratosphere  during  three  aircraft  campaigns  em¬ 
ploying  NASA’s  ER-2  aircraft.  The  three  campaigns 
are  (1)  the  second  Arctic  Airborne  Stratospheric  Ex¬ 
pedition  (AASE  II),  (2)  the  Stratospheric  Photochem¬ 
istry  Aerosol  and  Dynamics  Experiment  (SPADE),  and 
(3)  the  Airborne  Southern  Hemisphere  Ozone  Experi¬ 
ment/Measurements  for  Assessing  the  Effects  of  Strato- 


Table  1.  Aircraft  Missions 


Designation 

Dates 

Primary  Bases 

ER-2  Flights 

Stratosphere- Troposphere  Exchange  Project 
Tropical  Experiment  (STEP) 

Jan.-Feb.  1987 

Darwin,  Australia  (30°S,  131°E  ) 

17 

Airborne  Antarctic  Ozone  Experiment 
(AAOE 

Aug.-Oct.  1987 

Punta  Arenas,  Chile(53°S,  70°W  ) 

17 

Airborne  Aircraft  Stratospheric  Expedition 

(AASE  II) 

Dec.  1988  to 
March  1989 

Stavanger,  Norway(59°N,  6°E  ) 

16 

Airborne  Aircraft  Stratospheric  Expedition  II 
(AASE  II) 

Oct.  1991  to 
March  1992 

Fairbanks,  Alaska  (65°N,147°W) 

Bangor,  Maine  (45°N,  69°W) 

26 

Stratospheric  Photochemistry,  Aerosol  , 
and  Dynamics  Experiment  (SPADE) 

April-May  1993 

Moffet  Field,  California  (37°N,122°W) 

10 

Airborne  Southern  Hemisphere  Ozone 
Experiment  and  Measurements  for 

Assessing  the  Effects  of  Stratospheric 

Aircraft  (ASHOE/MAESA) 

Feb.-Nov.  1994 

Barbers  Point,  Hawaii  (21°N,158°W) 

Nadi,  Fiji  (18°S,177°E) 

Christchurch,  New  Zealand  (44°S,173°E) 

37 

References  describing  these  campaigns  are  STEP/Tropical  [Russell  et  al.,  1993],  AAOE  [Tuck  et  al.,  1989],  AASE 
[Turco  et  al.,  1990],  AASE  II  [Anderson  and  Toon,  1993] 


BACMEISTER  ET  AL.:  STRATOSPHERIC  HORIZONTAL  WAVENUMBER  SPECTRA 


9443 


spheric  Aircraft  (ASHOE/MAESA).  The  dates  and  lo¬ 
cations  of  these  and  three  earlier  aircraft  campaigns  em¬ 
ploying  the  ER-2  are  summarized  in  Table  1.  Maps 
showing  each  of  the  ER-2  flight  tracks  during  AASE  II, 
SPADE,  and  ASHOE/MAESA  are  displayed  in  Figure 
1. 

We  use  high  resolution  data  (sampling  frequency  1 
Hz)  from  the  approximately  70  ER-2  flights  conducted 
during  AASE  II,  SPADE,  and  ASHOE/MAESA.  A 
comparable  amount  of  data  was  available  from  the  three 
earlier  missions  (Table  1).  However,  improvements  in 
instrumentation  after  AASE  led  us  to  focus  on  the 
latest  three  missions.  To  our  knowledge,  this  is  the 
most  extensive  single  data  set  available  for  studies  of 
mesoscale  motions  in  the  middle  stratosphere.  The  ER- 
2’s  cruising  altitude  is  close  to  20  km  or  60,000  ft,  which 
is  20,000-30,000  ft  higher  than  the  commercial  airlin¬ 
ers  used  in  the  GASP  investigations.  Pressures  below 
70  mbar  are  reached  routinely,  and  measured  Brunt- 
Vaisala  frequencies  routinely  exceed  0.020  rad  s-1  [e.g. 
Gary ,  1989].  The  data  set  has  a  definite  bias  toward 
wintertime  polar  and  mid-latitude  conditions,  since  the 
primary  objective  of  four  of  the  six  missions  was  to  ex¬ 
amine  polar  ozone  chemistry.  However,  during  AASE 
II  and  SPADE  equatorward  flights  were  conducted  reg¬ 
ularly,  and  during  ASHOE/MAESA  two  out  of  every 
six  deployments  took  place  from  Hawaii  (see  Figure  1). 
Thus  approximately  one  third  of  the  flights  examined 
here  sampled  tropical  air. 

Velocity,  temperature,  and  pressure  data  were  mea¬ 
sured  from  the  ER-2  by  the  Meteorological  Measure¬ 
ment  System  (MMS)  [Scott  et  at .,  1990].  MMS  makes 
in  situ  measurements  of  the  three-dimensional  vector 
wind  velocity,  as  well  as  temperature  and  pressure,  at  a 
sampling  frequency  of  5  Hz.  Ozone  measurements  were 
obtained  by  the  NOAA  dual-beam  ozone  photometer 
[Proffit  et  al.,  1989]  at  a  sampling  frequency  of  1  Hz. 
We  also  examine  measurements  of  N2O  taken  by  the  air¬ 
borne  tunable  laser  absorption  spectrometer  (ATLAS) 
[Loewenstein  et  al.,  1989;  Podolske  et  al.,  1989].  AT¬ 
LAS  measurements  were  recorded  at  1  Hz.  However, 


AASE  — II  &  SPADE  ASHOE/MAESA 


Figure  1.  (Left)  Map  showing  ER-2  flight  tracks 
for  AASE  II  and  SPADE  aircraft  campaigns.  Flights 
during  these  two  campaigns  were  launched  from  Fair¬ 
banks,  Alaska;  Moffet  Field,  California;  and  Bangor, 
Maine.  (Right)  Map  showing  ER-2  flight  tracks  for 
ASHOE/MAESA  campaign.  Flights  originated  from 
Moffet  Field,  California;  Barbers  Point,  Hawaii;  Nadi, 
Fiji;  and  Christchurch,  New  Zealand. 


instrument  noise  appears  to  dominate  the  signal  for  fre¬ 
quencies  higher  than  0.3-0. 5  Hz  [Strahan  and  Mahlman, 
1994].  We  use  the  data  as  reported,  but  caution  is  re¬ 
quired  when  interpreting  the  highest  frequency  signals 
in  the  N2O  data 

For  purposes  of  characterizing  the  meteorological  con¬ 
ditions  in  which  measurements  were  taken,  we  used 
temperature  profile  data  from  the  microwave  temper¬ 
ature  profiler  (MTP)  [Gary,  1989;  Denning  et  al.,  1989] 
as  well  as  global  meteorological  analyses  of  wind  and 
temperature  obtained  from  the  National  Meteorological 
Center  Climate  Analysis  Center  (NMC/CAC)  [ Newman 
et  al,  1993].  The  MTP  provides  15  independent  tem¬ 
perature  measurements  within  a  layer  extending  about 
3  km  above  and  below  the  aircraft’s  altitude.  The  tem¬ 
perature  profiles  are  converted  to  potential  temperature 
using  pressures  derived  from  the  hydrostatic  relation 
[Gary,  1989].  A  complete  profile  is  obtained  every  14  s. 
The  NMC  data  are  obtained  daily  at  1200  UT  on  a 
global  grid  at  a  resolution  of  5°  x  2°  on  18  standard 
pressure  levels.  These  data  are  interpolated  in  all  three 
spatial  dimensions  as  well  as  in  time  onto  each  of  the 
ER-2  flight  tracks  flown  during  AASE  II,  SPADE,  and 
ASHOE/MAESA.  Figure  1  shows  maps  of  each  of  the 
ER-2  flight  tracks  executed  during  AASE  II,  SPADE, 
and  ASHOE/MAESA. 

2.2  Wavelet  Transforms 

Our  analysis  of  the  aircraft  data  employs  wavelet 
transforms  rather  than  Fourier  transforms.  Orthonor¬ 
mal  wavelets  are  families  of  orthogonal  functions  which, 
unlike  sines  and  cosines,  are  localized  in  both  the  spa¬ 
tial  and  wavenumber  domains  [Daubechies,  1988;  Mal- 
lat,  1989].  A  large  number  of  families  of  orthonormal 
wavelets  has  been  identified.  Some  are  compactly  sup¬ 
ported  (i.e.  exactly  equal  to  zero  outside  of  a  finite  in¬ 
terval),  while  others  approach  zero  asymptotically  out¬ 
side  of  a  finite  interval.  In  the  following  analysis,  we 
use  the  discreet  wavelet  transform  (DWT)  described  by 
Press  et  al.  [1992].  This  algorithm  can  accommodate 
any  compactly  supported  wavelets  for  which  filter  coef¬ 
ficients  are  known.  We  use  the  six-coefficient  member 
of  the  wavelet  class  identified  by  Daubechies  [1988].  An 
example  of  this  wavelet  can  be  seen  in  Figure  2a. 

Due  to  the  localization  of  the  basis  functions,  DWT 
spectra  contain  information  about  the  location  of  vari¬ 
ability  in  the  time  series  as  well  as  its  frequency.  The 
DWT  spectrum  is  most  easily  interpreted  when  bro¬ 
ken  down  into  its  component  “octaves” .  At  each  octave 
the  DWT  shows  variability  at  a  specific  frequency  or 
wavenumber  as  a  function  of  location.  Adjacent  octaves 
show  frequencies  separated  by  factors  of  2.  Figures  2b- 
2d  show  the  DWT  of  the  “data”  in  Figure  2a  broken 
down  into  some  of  its  component  octaves.  An  obvious 
relationship  exists  between  the  peaks  in  the  DWT  spec¬ 
trum  and  the  localized  variability  in  the  data.  Unfor¬ 
tunately,  for  a  data  sequence  with  2N  elements,  only  N 
octaves  can  be  obtained.  Thus,  the  price  for  the  spatial 
information  contained  in  the  DWT  is  a  significant  loss 
of  frequency  resolution  when  compared  to  a  fast  Fourier 
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Original  data 


X. 


i/=9  (2"=512) 


x. 


Figure  2.  (a)  Artificial  data  sequence  consisting  of  a  high-frequency  Gaussian-modulated  pulse 

described  by  g-C*-012)2/0  022[sin(2007ra:)-|-sin(10007ra:)]  and  a  Daubechies  (6)  wavelet  at  x  =0.35- 
0.5.  The  wavelet  is  the  inverse  transform  of  the  unit  vector  e 44  ,  i.e. ,  zeros  everywhere  except  for 
the  44th  element,  which  is  1.  Since  44  lies  between  25  =  32  and  26  =  64  the  peak  in  the  forward 
transform  occurs  in  fifth  octave  at  a  position  close  to  x  =  3/8  ±  3/32,  where  the  uncertainty 
results  from  arbitrariness  in  the  positioning  of  the  six  element  filter,  (b)  The  fifth  octave  (v  =  5) 
of  the  DWT  of  the  data  in  Figure  2a.  (c)  As  in  Figure  2b,  except  for  v  =  7.  (d)  As  in  Figure  2b, 
except  for  v  =  9. 


transform(FFT).  However,  when  plotted  on  logarithmic 
axes,  as  we  will  do  later,  the  DWT  frequencies  are  eq- 
uispaced,  and  the  degraded  frequency  resolution  is  no 
longer  a  serious  drawback. 

The  advantages  of  DWTs  over  FFTs  are  greatest  for 
sequences  of  data  which  are  characterized  by  large  inter¬ 
mittent  events.  In  this  case,  DWTs  provide  a  straight¬ 
forward  connection  between  features  in  the  original  data 
and  spectral  features.  This  is  clearly  seen  in  Figure 
2.  DWTs  were  chosen  for  this  analysis  because  we  ex¬ 
pected  the  aircraft  data  to  consist  primarily  of  inter¬ 
mittent  high-energy  peaks  in  variance  associated  with 
isolated  wave  or  turbulence  events.  This  supposition  is 
based  on  previous  studies  of  atmospheric  wind  measure¬ 
ments  using  aircraft  [e.g .,Nastrom  et  al.,  1987;  Yamada 


and  Ohkitani,  1991;  Nastrom  and  Fritts,  1992;  Bacmeis- 
ter  et  al.,  1994]. 

DWTs  also  offer  several  procedural  advantages  over 
FFT’s.  Application  of  windows  to  the  data  to  mitigate 
edge  effects  is  unnecessary.  Edge  effects  can  be  sim¬ 
ply  accounted  for  by  ignoring  DWT  components  located 
near  the  edge  of  the  original  data  sequence.  Similarly, 
the  possibility  of  artificial  variability  induced  by  aircraft 
maneuvers  (e.g.,  dives,  box  turns)  is  easily  assessed  by 
examination  of  the  DWT  in  the  vicinity  of  the  maneu¬ 
vers.  However,  it  should  be  emphasized  that  none  of 
the  results  presented  in  this  study  depend  crucially  on 
the  use  of  DWTs  instead  of  FFTs,  nor  does  the  partic¬ 
ular  choice  of  wavelet  function  significantly  affect  the 
results. 
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2.3  Determination  of  Power  Spectral  Densities 
From  DWTs 

Results  of  earlier  studies  of  small-scale  variability  in 
the  atmosphere  are  typically  presented  in  terms  of  av¬ 
eraged  power  spectra.  The  most  important  quantity 
derived  from  these  spectra  is  usually  the  mean  slope  of 
the  spectra  in  log-log  coordinates,  i.e.,  the  spectral  in¬ 
dex.  For  comparison  of  our  DWT  spectra  with  Fourier 
power  spectra  presented  in  earlier  studies,  we  must  first 
define  a  mean  DWT  power  spectrum  for  our  aircraft 
data  as  a  function  of  frequency  or  wavenumber.  We  do 
this  following  the  procedure  outlined  by  Yamada  and 
Ohkitani  [1991]  and  Katul  et  al.  [1994]. 

A  typical  ER-2  flight  during  the  AASE  II,  SPADE, 
or  ASHOE/MAESA  missions  was  6-8  hours  long.  The 
input  data  for  the  DWT  analysis  consist  of  5-Hz  MMS 
temperature,  pressure  and  wind  velocity  measurements, 
which  have  been  bin-averaged  onto  a  uniform  1-s  time 
grid,  and  the  original  1-Hz  NOAA  ozone  and  ATLAS 
N2O  measurements.  Thus  the  complete  time  series  of 
each  quantity  for  each  flight  is  typically  ~25,000  ele¬ 
ments  long.  Infrequent  data  gaps  were  filled  using  linear 
interpolation. 

We  calculate  average  DWT  power  spectra  for  subin¬ 
tervals  of  data  lasting  1024  s,  so  that  20  to  25  spectra 
are  typically  obtained  for  each  quantity  for  each  flight. 
First,  the  DWT  spectrum  is  calculated  for  the  jth  1024- 
element  data  segment,  Xj(t), 

Xj  ( t )  — ►  DWT  — ►  [Djt 0,  Dh  1, ...,  Dj: 9]  =  [Dyv], 
where  1/  —  0,...,9.  The  Dj^u  denote  the  10  separate 


octaves  v  of  the  DWT  of  Xj(t),  each  containing  2"  el¬ 
ements.  These  are  used  to  calculate  an  average  DWT 
squared-amplitude  within  each  octave  for  Xj(t), 

WTiM=(^)ZDU‘)-  (i) 

'  '  t=i 

A  mean-squared  wavelet  amplitude  is  thus  defined  for 
each  of  the  10  octaves  which  can  be  obtained  from 
a  1024-element  time  series.  The  index  of  the  octave, 
u,  corresponds  roughly  to  log2(Q),  where  D  is  the  fre¬ 
quency  of  the  FFT  spectrum.  Fourier  modes  with  fre¬ 
quencies  D  =  2"  have  DWT  power  spectra  which  peak 
sharply  at  the  fcth  octave.  Fourier  modes  for  which 
2"-1  <  Q  <  2"  project  power  onto  the  two  bracketing 
octaves  u  and  v  —  1. 

The  total  number  Nap  of  individual  DWT  power  spec¬ 
tra  obtained  for  each  measured  quantity  varies  between 
1290  and  1680,  when  the  entire  data  set  of  73  or  so 
flights  is  used.  Differences  in  the  number  of  individual 
spectra  for  the  various  quantities  are  due  to  occasional 
payload  changes,  varying  on-off  times  for  different  in¬ 
struments,  data  reduction  failures,  and  instrument  fail¬ 
ures.  The  more  than  1000  individual  DWT  spectra  for 
each  quantity  can  be  log-averaged  to  give  an  average 
power  spectrum  for  the  entire  data  set, 

log ,[vrr*V)]  =  (^)  ]r>g2[ wts[v)].  (2> 

This  average  DWT  power  spectrum  is  comparable  to 
the  average  FFT  power  spectral  densities  presented  in 


Power  Law  Time  Series 


Comparison  of  DWT  and  FFT 


Figure  3.  (top)  Segment  of  an  artificial  power  law  time  series  with  energy  or  variance  3, 
(bottom)  Comparison  of  power  spectral  densities  (PSD)  for  the  data  in  Figure  3  (top)  obtained 
from  a  windowed  FFT  (solid  line)  and  a  DWT  (solid  circles).  The  expected  -3  slope  is  indicated 
by  the  thick  dashed  line. 
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other  studies.  Average  power  spectra  for  subsets  of  the 
data  are  calculated  in  an  analogous  manner. 

In  order  to  test  this  procedure,  we  applied  it  to  a 
1024-element  subinterval  of  an  artificially  generated, 
4096-element  aircraft  time  series, 

2048 

Y(t)=  Y  Q~15  sin[Qt  +  <ft(Q)] 

n=i 

where  t  represents  the  interval  [0,2tt]  divided  into  4096 
segments  and  <f>(Ll)  is  a  random  phase  between  0  and  2 it. 
Results  are  shown  in  Figure  3.  An  FFT  power  spectrum 
of  the  same  time  series  is  shown  for  comparison.  The 
spectral  index,  or  slope  of  the  DWT  power  spectrum 
in  log-log  coordinates,  agrees  well  with  the  analytical 
spectral  index  of  —3.  The  average  slope  of  the  FFT 
power  spectrum  in  log-log  coordinates  is  also  close  to 
—3.  The  scatter  or  “structure”  visible  in  the  FFT  spec¬ 
trum  results  entirely  from  the  application  of  a  triangu¬ 
lar  (Bartlett)  window  to  the  data  before  calculating  the 
power  spectrum.  This  is  necessary  in  practice  to  reduce 
contamination  of  the  computed  FFT  spectrum  by  edge 
effects. 

In  the  following  analysis,  the  approximate  mean  ER-2 
airspeed  of  200  ms-1  is  used  to  transform  timescales  in 
the  data  into  spatial  scales  via  a  straightforward  Tay¬ 
lor  transformation.  While  a  standard  method  in  the 
presence  of  turbulence,  Taylor’s  hypothesis  is  not  auto¬ 
matically  valid  in  the  presence  of  wave  motions,  which 
propagate  relative  to  the  background  wind.  In  section 
4  it  will  be  shown  that  Taylor’s  hypothesis  does  not  in¬ 
troduce  significant  errors  for  the  gravity  wave  spectra 
that  are  typically  encountered  at  these  heights.  Thus 
we  assume  that  a  feature  in  the  aircraft  time  series  with 
a  period  of  1024  s  possesses  a  horizontal  wavelength  of 
204.8  km,  while  a  feature  at  our  resolution  limit  of  2 
s  possesses  a  wavelength  of  400  m.  For  direct  compar¬ 
ison  with  FFT  power  spectral  densities  (PSDs)  shown 
in  previous  studies,  we  scale  the  mean-squared  wavelet 
amplitudes  in  (1)  and  (2)  as  in  the  work  by  Katul  et  al. 
[1994]: 

where  WT{u)  can  represent  either  the  individual  DWT 
squared-amplitudes  in  (1)  or  the  log  -  averaged  ampli¬ 
tudes  in  (2)  and  As  represents  the  distance  covered  by 
the  aircraft  between  samples.  We  use  As  »  200  m. 
E{v)  is  now  directly  comparable  to  Fourier  power  spec¬ 
tral  densities  calculated  in  earlier  studies  of  atmospheric 
horizontal  wavenumber  spectra. 

3.  Results 

3.1  Average  Spectra  of  Velocities  and  Potential 
Temperature 

Figure  4a  shows  DWT  power  spectral  densities  E(v) 
of  horizontal  and  vertical  velocities,  log-averaged  over 
the  entire  data  set.  The  set-averaged  power  spectrum 


of  potential  temperature  is  shown  in  Figure  4b.  Spec¬ 
tral  indices  (slopes)  are  significantly  steeper  than  —5/3 
for  scales  shorter  than  6  km  ( v  >  5)  for  all  wind  com¬ 
ponents  as  well  as  for  potential  temperature.  The  spec¬ 
tral  slopes  at  these  scales  appear  to  be  closer  to  —2.5  or 
—3.  At  wavelengths  between  6  and  25  km  the  spectrum 
of  both  horizontal  velocity  components  flattens,  and  at 
larger  scales  approaches  the  more  familiar  —5/3  law  re¬ 
ported  in  studies  using  the  GASP  data  set  [Nastrom  et 
al.,  1984;  Nastrom  and  Gage,  1985]  and  other  studies 
using  radar  [e.g.  Gage  and  Nastrom,  1985].  Analyses 
of  higher-resolution  GASP  data  [ Nastrom  et  al. ,  1984; 
Nastrom  and  Gage,  1985]  returned  indices  of  approxi¬ 
mately  —5/3  down  to  scales  of  2  km. 

Thus  the  average  spectra  in  Figure  4  appear  to  dis¬ 
agree  with  the  GASP  results  at  horizontal  scales  be¬ 
low  about  25  km.  Later,  more  detailed  analysis  by 
Nastrom  et  al.  [1987]  found  increased  spectral  pow¬ 
ers  and  spectral  indices  near  -3  at  large  wavenumbers 
in  GASP  spectra  measured  over  mountainous  terrain  in 
high  wind  conditions.  The  GASP  spectra  under  these 
conditions  are  similar  to  the  set-averaged  spectra  shown 
in  Figure  4.  However,  the  GASP  climatological  mean 
spectrum  under  normal  meteorological  conditions  re¬ 
tained  a  spectral  index  close  to  -5/3  down  to  scales 
of  2  km  [Nastrom  et  al.,  1987],  The  average  spectra 
shown  in  Figure  4  contain  measurements  made  in  a  va¬ 
riety  of  meteorological  conditions.  As  can  be  seen  in 
Figure  1,  a  large  fraction  (in  fact  well  over  half)  of  the 
ER-2  flight  time  for  which  we  have  data  occurred  over 
oceans.  Thus  it  appears  unlikely  that  the  agreement 
between  our  average  spectra  and  the  high-wind,  high- 
topography  spectra  examined  by  Nastrom  et  al.  [1987] 
is  due  to  a  topographic  bias  in  our  data.  In  section 
3.3  we  will  examine  in  more  detail  the  effect  of  several 
meteorological  variables  on  our  average  spectra. 

At  high  wavenumbers  the  average  velocity  spectra  in 
Figure  4a  resemble  some  of  those  derived  from  other 
high-resolution  data  sets  obtained  using  specially  in¬ 
strumented  aircraft  [Shur,  1962,  1971;  Vinnichenko  et 
al.,  1965;  Koprov,  1965;  Myrup,  1968;  Lilly  and  Lester, 
1974].  Vinnichenko  et  al.  [1965]  and  Shur  [1971]  de¬ 
scribed  a  series  of  experiments  conducted  over  central 
Asia  (Uzbekistan)  [ Vinnichenko  et  al.,  1965]  and  the 
Caucasus  Mountains  [Shur,  1971].  In  these  studies  air¬ 
craft  were  flown  at  several  altitudes  between  the  lower 
troposphere  and  the  tropopause.  Horizontal  velocity 
spectra  with  spectral  indices  between  —5/3  and  —3  were 
obtained  for  scales  down  to  500  m.  The  steeper  spec¬ 
tral  slopes  were  invariably  associated  with  high  ambient 
stratification.  Lilly  and  Lester  [1974]  described  a  series 
of  flights  over  the  Sangre  de  Cristo  Mountains  in  Col¬ 
orado  at  altitudes  of  13-20  km.  All  of  the  horizontal 
velocity  spectra  measured  during  these  flights  had  in¬ 
dices  of  approximately  —3  down  to  the  reported  resolu¬ 
tion  limit  of  2  km.  Koprov  [1965]  and  Myrup  [1968] 
described  separate  aircraft  experiments  conducted  in 
strong,  low-level,  nocturnal  inversions  over  flat  desert 
terrain  in  central  Asia  [Koprov,  1965]  and  over  a  dry 
lake  bed  in  northern  California  [Myrup,  1968].  In  each 
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Figure  4.  Log-averaged  DWT  power  spectra  of  (a)  velocities  and  (b)  potential  temperature. 
Log  averages  are  taken  over  the  entire  set  of  available  1024-s  spectra  of  each  quantity,  (a)  Solid 
circles  show  log-averaged  spectrum  of  zonal  velocity,  solid  diamonds  show  meridional  velocity, 
and  solid  triangles  show  vertical  velocity.  Short  dashed  shows  a  —3  slope.  Long  dashed  line 
shows  a  —5/3  slope. 


of  these  experiments,  kinetic  energy  spectra  with  in¬ 
dices  near  —3  were  measured  during  the  early  morn¬ 
ing  while  the  inversion  was  strong.  As  the  inversion 
was  destroyed  by  solar  heating,  the  indices  of  the  mea¬ 
sured  kinetic  energy  spectra  approached  —5/3.  Finally, 
Berman  [1976]  reported  horizontal-velocity  frequency 
spectra  ~fi-3  measured  by  a  ship-based,  tethered  bal¬ 
loon  at  the  top  of  a  stable  marine  boundary  layer. 

In  each  of  the  studies  just  described,  spectra  with 
slopes  steeper  than  -5/3  were  obtained  only  in  atypical 
geographical  or  meteorological  environments.  However, 
the  stratospheric  meteorological  data  collected  from  the 


ER-2  appear  to  yield  spectra  with  slopes  close  to  —2.5 
under  general  conditions.  A  randomly-chosen  subset  of 
individual  power  spectra  for  1024-s  segments  of  zonal 
velocity  (17)  data  is  shown  in  Figure  5.  The  individual 
spectra  tend  to  lie  within  a  fairly  narrow  band  around 
the  average  power  spectrum  of  U  in  Figure  4a.  To  make 
this  characterization  more  quantitative  and  to  assess 
the  statistical  significance  of  our  slope  estimates,  we  ex¬ 
amine  the  distribution  of  slope  values  obtained  at  each 
octave  for  each  of  the  N,p  >  1000  individual  1024-s 
DWT  power  spectra  of  each  quantity. 

The  value  of  the  spectral  index  between  the  vth.  oc- 
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Figure  5.  Randomly  selected  subset  of  individual 
DWT  power  spectra  for  1024-s  segments  of  zonal  wind 
data  from  MMS.  Note  the  relatively  tight  concentration 
of  the  lines.  This  suggests  that  the  mean  spectrum  in 
Figure  4a  is  robust. 


tave  and  the  (v  +  l)th  octave  for  each  individual  spec¬ 
trum  is  given  by, 

Sj  {v  +  1/2)  =  log \2[Ej(v  +  1)]  -  \og2\Ej(y)}  (3) 

where  Ej{v)  is  the  average  power  of  the  uth  octave  of 


the  DWT  for  the  jth  data  segment  as  defined  in  (1). 
There  is  a  denominator  implicit  in  (3)  which  measures 
the  change  in  log-frequency  or  log- wavenumber,  Q,  be¬ 
tween  octaves  u+ 1  and  i/,  but  is  simply  equal  to  1  since 
D(i/)  ~  2".  Figure  6  shows  histograms  of  slope  values 
Sj  of  U  at  each  octave  for  the  entire  data  set.  The 
histograms  are  sharply  peaked,  particularly  at  higher 
octaves  .  The  histograms  for  u  +  1/2  =  5.5, 6.5  and 
7.5  have  peaks  between  —2.33  and  —2.66  with  few  oc¬ 
currences  of  slopes  flatter  than  —5/3.  The  histograms 
strongly  suggest  that  the  mean  spectral  slopes  for  hori¬ 
zontal  velocity  in  this  data  set  at  scales  below  25  km  are 
in  fact  significantly  steeper  than  —5/3.  Table  2  shows 
mean  values  and  80%  confidence  limits  for  the  slope  his¬ 
tograms  of  U  in  Figure  6,  as  well  as  confidence  limits 
for  slope  histograms  of  the  other  quantities  examined. 

The  mean  power  spectrum  of  vertical  velocity  (Figure 
4a)  differs  markedly  from  those  of  the  other  meteoro¬ 
logical  quantities  examined  here.  At  scales  of  3  km  and 
shorter  (v  >  6)  the  slope  of  the  spectrum  is  near  —2.4, 
which  is  similar  to  the  slope  of  the  horizontal  veloc¬ 
ity  spectra  at  these  scales.  This  agreement  would  be 
consistent  with  an  explanation  based  on  homogeneous 
buoyancy-modified  turbulence  in  which  the  power  spec¬ 
trum  of  all  velocity  components  decays  as  k~r ,  where 
r  is  between  4/3  and  3,  [e.g.,  Lumley,  1964;  Weinstock, 
1978,  1985a].  In  fact,  vertical  velocity  spectra  with  ob¬ 
served  k~3  behavior  at  scales  0.5-5  km  were  reported 
by  Koprov  [1965],  Myrup  [1968],  and  Shur  [1971]. 

However,  at  scales  longer  than  6  km  the  slope  of  our 
average  power  spectrum  for  vertical  velocity  flattens  out 
rapidly  and  above  25  km  the  spectrum  is  essentially 
flat.  The  histogram  analysis  of  the  slopes  (Table  2) 
again  suggests  that  slope  estimates  are  relatively  ro¬ 
bust.  This  flattening  out  of  the  vertical  velocity  power 
spectrum  at  large  scales  simply  reflects  the  intuitively 
obvious  fact  that  strong,  mesoscale  updrafts  and  down- 
drafts  do  not  exist  in  the  mean  atmosphere.  However, 
to  our  knowledge  this  necessary  feature  of  the  vertical 
velocity  wavenumber  spectrum  has  not  been  observed 
by  previous  aircraft  studies. 

Interestingly,  ST  radar  measurements  of  vertical  ve¬ 
locities  yield  spectra  with  spectral  indices  close  to  0  at 
low  frequencies  [Ecklund  et  al.,  1986;  VanZandt  et  al, 
1991].  In  these  studies,  the  vertical  velocity  spectrum 
typically  posesses  a  steep  high-frequency  tail,  while  ex¬ 
hibiting  nearly  constant  power  (i.e. ,  spectral  index  0  ) 
for  measured  periods  longer  than  10-30  min.  A  rigorous 
treatment  of  these  spectra  involves  the  consideration 
of  Doppler  shifting  effects  and  Taylor  transformations, 
in  order  to  identify  their  dynamical  content  [ Gage  and 


Figure  6.  Histograms  showing  distribution  of  slopes  for  each  octave  for  all  individual  horizontal 
wind  DWT  power  spectra  in  the  data  set.  (a)  Histogram  of  slope  values  between  octaves  1  and 
2,  or  roughly  speaking,  for  horizontal  scales  between  102.4  km  and  51.2  km.  (b)  Same  as  Figure 
6a  for  octaves  2  and  3,  and  horizontal  scale  between  51.2  and  25.6  km  .  (c)  Octaves  3  and  4, 
horizontal  scale  between  25.6  and  12.8  km.  (d)  Octaves  4  and  5,  horizontal  scale  between  12.8 
and  6.4  km.  (e)  Octaves  5  and  6,  horizontal  scale  between  6.4  and  3.2  km.  (f)  Octaves  6  and  7, 
horizontal  scale  between  3.2  and  1.6  km.  (g)  Octaves  7  and  8,  horizontal  scale  between  1.6  km 
and  800  m.  (h)  Octaves  8  and  9,  horizontal  scale  between  800  m  and  400  m. 
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Table  2.  Spectral  Indices 


Octave 

10th  Percentile 

Mean 

90th  Percentile 

n 

2.5 

+0.17 

Zonal  Velocity 
-1.66 

-3.17 

1348 

4.5 

-0.83 

-1.99 

-2.83 

6.5 

-1.83 

-2.51 

-2.83 

2.5 

+0.50 

Meridional  Velocity 
-1.70 

-3.17 

1348 

4.5 

-0.83 

-1.99 

-2.83 

6.5 

-1.83 

-2.55 

-2.83 

2.5 

+1.83 

Vertical  Velocity 
-0.02 

-3.33 

1290 

4.5 

+0.17 

-0.94 

-1.83 

6.5 

-1.50 

-2.43 

-2.83 

2.5 

+0.17 

Potential  Temperature 
-1.91 

-3.50 

1510 

4.5 

-0.83 

-2.07 

-2.83 

6.5 

-1.83 

-2.52 

-2.83 

2.5 

+0.50 

Ozone 

-1.78 

-3.50 

1639 

4.5 

-0.50 

-1.74 

-2.83 

6.5 

-1.17 

-2.04 

-2.83 

2.5 

+0.17 

N20 

-1.81 

-3.50 

1160 

4.5 

-0.50 

-1.74 

-2.83 

6.5 

-0.83 

-1.90 

-2.83 

Results  of  histogram  analysis  of  spectral  slopes  for  each  quantity  examined.  Slope  histograms  for  horizontal 
velocity  are  shown  in  Figure  5.  The  table  shows  the  10th  and  90th  percentile  values  for  the  slope  histograms  as  well 
as  the  mean  value  for  v  =  2.5,  4.5,  and  6.5.  The  last  column  gives  the  total  number  of  1024-s  data  segments  available 
for  each  quantity. 


Nastrom,  1985;  Vincent  and  Eckermann,  1990;  Van- 
Zandt  et  al.,  1991].  A  similar  sort  of  analysis  of  the 
aircraft  spectra  is  undertaken  in  Section  4. 

3.2  Average  Spectra  of  Ozone  and  N2O 

Figure  7  shows  set-averaged  spectra  of  ozone  mixing- 
ratio  variance  and  N2O  mixing-ratio  variance.  They  are 
noteworthy  in  that  they  are  the  only  spectra  analyzed 
which  appear  to  conform  closely  to  the  —5/3  power  law 
at  most  scales.  Slope  histograms  yield  spectral  indices 
close  to  —1.7  (Table  2).  A  relatively  distinct  transi¬ 
tion  from  a  —5/3  behavior  to  —2  behavior  occurs  in 
the  ozone  spectrum  at  a  wavelength  of  3.2  km  (u  =  6). 
However,  the  slope  histograms  (not  shown)  suggest  that 
this  is  due  more  to  a  reduction  in  the  relative  number 
of  spectral  segments  with  flat  slopes,  rather  than  to  a 
shift  in  the  peak  of  the  distribution  of  slopes.  The  peaks 
remain  near  —5/3  for  all  octaves,  except  for  the  highest 


(1/  +  1/2  =  8.5)  octave,  which  may  be  an  aliasing  ef¬ 
fect  [Eckermann,  1990].  It  is  also  worth  noting  that  the 
spread  in  the  slope  histograms  for  ozone  is  somewhat 
larger  than  for  horizontal  velocity  and  vertical  velocity. 
This  is  seen  in  the  slightly  larger  spread  between  the  10 
and  90  percentiles  of  the  spectral  index  determinations 
for  ozone  and  horizontal  velocity  given  in  Table  2. 

The  average  power  spectrum  for  N2O  conforms  closely 
to  k~ 5/,s  at  all  scales.  This  is  consistent  with  the  anal¬ 
ysis  of  Strahan  and  Mahlman  [1994],  who  derived  spec¬ 
tral  indices  between  -2  and  -5/3  for  N2O  data  from  the 
AASE  and  AAOE  campaigns.  Thus  the  ozone  spectra 
and  the  N2O  spectra  derived  from  the  present  data  ex¬ 
hibit  similar  k~p  behavior  where  1.5  <  p  <  2.0  for  the 
range  of  horizontal  scales  100  km  >  A  >  1  km.  Thus 
the  average  power  spectra  of  tracer  variance  appear  to 
be  significantly  flatter  than  the  power  spectra  of  the 
dynamical  quantities  at  these  scales. 


BACMEISTER  ET  AL.:  STRATOSPHERIC  HORIZONTAL  WAVENUMBER  SPECTRA 


9451 


wavenumber(cyc  m~’) 
1CT5  lO-*  I  O’’ 


0  123456789  10 

Octave  (i/) 

Figure  7.  Log-averaged  DWT  power  spectra  of  ozone 
and  N2O  mixing  ratio  variance.  Log  averages  are  taken 
over  the  entire  set  of  available  1024  s  spectra  of  each 
quantity.  Solid  circles  show  log-averaged  spectrum  of 
ozone  mixing  ratio  variance.  Solid  triangles  show  N2O 
variance.  Mixing  ratios  N2O  were  arbitrarily  multiplied 
by  8  to  make  them  comparable  to  ozone  mixing  ratios. 
Short  dashed  line  shows  —3  slope.  Long  dashed  line 
shows  —5/3  slope. 


3.3  Spectra  Stratified  According  to  Velocity  Vari¬ 
ance,  Mean  Velocity,  and  Stratification 

Several  mechanisms  can  be  invoked  to  explain  the 
steeper  than  expected  slopes  of  our  set-average  dynam¬ 
ical  power  spectra.  The  possibility  of  turbulence  in  the 
buoyant  subrange  has  already  been  mentioned.  Other 
possibilities  include  spectra  due  to  saturated  or  unsat¬ 


urated  gravity  wave  motions.  In  order  to  obtain  possi¬ 
ble  insights  into  the  nature  of  the  motions  creating  our 
power  spectra,  we  now  examine  the  impact  of  several 
meteorological  variables  on  the  shape  of  the  spectra  ob¬ 
tained  from  our  aircraft  data.  A  similar  classification 
of  GASP  spectra  according  to  mean  ambient  horizontal 
velocity  was  performed  by  Nastrom  et  al.  [1987]. 

Here  we  examine  not  only  the  influence  of  the  mean 
ambient  horizontal  velocity  on  the  shape  of  the  mea¬ 
sured  spectra,  but  also  that  of  mean  horizontal  velocity 
variance  and  mean  stratification.  Mean  stratification 
is  expected  to  affect  the  spectral  slope  if  buoyancy- 
dominated  turbulence  is  present  [  Weinstock ,  1978;  Vin¬ 
nichenko  et  al.,  1980].  Mean  horizontal  velocity  and 
mean  horizontal  velocity  variance  may  be  expected  to 
influence  the  shape  of  a  spectrum  of  interacting  gravity 
waves  [Fritts  and  VanZandt,  1987;  Hines,  1991;  Ecker- 
mann,  1995],  as  well  as  a  turbulent  spectrum  in  the  2-D 
or  buoyant  subranges  [Weinstock,  1978;  Nastrom  et  al., 
1987;  Eckermann,  1990]. 

Estimates  of  mean  horizontal  velocity  and  rms  hori¬ 
zontal  velocity  variance  are  obtained  directly  from  the 
MMS  data.  Mean  ambient  horizontal  velocity  is  calcu¬ 
lated  for  each  1024-s  interval  in  the  data  set  by  simply 
bin  averaging  the  MMS  measured  horizontal  velocity 
over  the  interval. 

An  index  of  rms  horizontal  velocity  variance  is  ob¬ 
tained  directly  from  the  DWT  power  spectra  of  hori¬ 
zontal  velocity.  We  simply  take  the  log-average  power 
over  all  nine  octaves  of  the  DWT  power  spectrum  of 
horizontal  velocity  for  a  given  1024-s  interval  of  data. 
This  average  is  compared  to  the  average  over  the  entire 
data  set  to  characterize  unusually  active  or  quiet  peri¬ 
ods.  Estimates  of  the  mean  stratification  at  the  time  of 
measurement  are  taken  from  NMC  analyses  of  poten¬ 
tial  temperature  interpolated  to  the  location  and  time 
of  the  measurement  or  from  time-averaged  MTP  mea¬ 
surements  of  lapse  rates. 

Effect  of  horizontal  velocity  variance  To  sort 
spectra  according  to  total  horizontal  velocity  variance, 
first  we  calculate  the  log-averaged  PSD  for  total  hor¬ 
izontal  velocity  over  all  octaves  for  the  set-averaged 
power  spectrum  shown  in  Figure  4a. 

iog2  (rev)  =  Q)  Eiog2 

'  '  v  =  l 

where  the  set-averaged  spectrum  of  total  horizontal 
velocity  is  given  by  the  sum  of  the  components,  i.e., 
Euv  =  Ev  +  Ev  .  Set  averages  are  defined  by 
(2).  From  the  log-averaged  power  spectra  of  zonal 
and  meridional  velocity  shown  in  Figure  4a,  we  obtain 
»2x  25'6  m3  rad-1  s-2.  We  now  calculate  an 
analogous  quantity  \Ej\uv  for  each  of  the  Nsp  >  1300 
individual  DWT  power  spectra  of  horizontal  velocity. 
We  divide  the  data  into  two  groups:  an  “active”  group, 
in  which  [Ej]uv  >  2 [E  ]uv,  and  a  “quiet”  group,  in 
which  [Ej]uv  <  \[E  ]uv- 
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Figure  8.  DWT  power  spectra  of  (a)  horizontal  and  (b)  vertical  velocity  sorted  according  to  the 
magnitude  of  the  horizontal  velocity  variance  (see  text).  Solid  triangles  show  average  for  high 
variability  group.  Solid  circles  show  low  variability  group. 


Figure  8  shows  DWT  power  spectra  of  zonal  and  ver¬ 
tical  velocity  for  both  groups.  Power  spectra  of  verti¬ 
cal  velocity  are  shown  in  Figure  8b.  The  upper  curve 
in  Figure  8b  shows  the  log-averaged  power  spectrum 
of  vertical  velocity  for  active  intervals,  and  the  lower 
curve  shows  the  log-averaged  spectrum  for  quiet  inter¬ 
vals.  The  two  curves  are  simply  offset  in  the  vertical, 
and  no  detectable  difference  in  the  shape  of  the  verti¬ 
cal  velocity  spectra  appears  to  result  from  grouping  the 
data  in  this  way. 

Similar  sorting  of  the  horizontal  velocity  spectra  in 
Figure  8a  yields  a  similar  result.  A  mean  vertical  offset 
between  the  two  curves  in  Figure  8a  is  an  unavoidable 


result  of  the  sorting  procedure.  However,  there  is  no  re¬ 
quirement  that  spectral  slopes  remain  the  same.  Never¬ 
theless,  the  log-averaged  power  spectra  again  appear  to 
be  nearly  identical  in  shape  except  for  the  vertical  off¬ 
set.  This  result  appears  to  be  somewhat  at  odds  with 
the  theory  of  Weinstock  [1978],  which  suggests  that  the 
shape  of  the  kinetic  energy  spectrum  for  turbulence  in 
the  buoyant  subrange  should  depend  on  the  local  rms 
velocity  variance. 

Effect  of  mean  horizontal  velocity  The  situa¬ 
tion  changes  somewhat  when  the  spectra  are  sorted  ac¬ 
cording  to  the  mean  horizontal  wind  speed  at  the  time 
of  measurement.  Again,  we  sort  the  individual  spectra 
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Figure  9.  Log-averaged  DWT  power  spectra  of  (a)  horizontal  and  (b)  vertical  wind  sorted 
according  to  the  magnitude  of  mean  ambient  horizontal  velocity  (i.e.,  mean  horizontal  wind 
speed).  Mean  ambient  horizontal  velocity  is  obtained  from  1024  s  averages  of  ( U 2  +  V2)1/2  as 
measured  by  the  MMS  along  the  ER-2  flight  track.  Two  subsets  of  the  data  are  constructed:  one 
for  mean  horizontal  velocity  <10  m  s-1  and  another  for  mean  horizontal  velocity  >30  m  s-1. 
Solid  triangles  show  average  for  high  wind  speed  group.  Solid  circles  show  low  wind  speed  group. 


into  two  groups;  one  for  which  the  corresponding  1024- 
s-averaged  MMS  horizontal  wind  speed  is  greater  than 
30  m  s_1  and  another  for  which  it  is  less  than  10  ms'1. 
These  limits  were  chosen  to  yield  groups  of  spectra  with 
similar  numbers  of  members.  This  sorting  procedure  is 
similar  to  that  performed  by  Nastrom  et  al.  [1987]  with 
the  GASP  data,  except  that  they  used  a  single  velocity 
value  of  25  ms-1  to  separate  their  data  into  high-  and 
low-velocity  groups.  In  addition,  unlike  Nastrom  et  al. 
we  do  not  segregate  our  data  according  to  underlying 
topography. 


Figure  9b  shows  log-averaged  power  spectra  of  verti¬ 
cal  velocity  for  subsets  of  data  which  have  been  sorted 
according  to  mean  horizontal  velocity.  In  contrast  to 
the  earlier  subdivision  by  rms  horizontal  velocity  vari¬ 
ance,  the  shape  of  the  vertical  velocity  spectra  does  ap¬ 
pear  to  change  in  response  to  varying  mean  horizontal 
wind  speeds.  In  high  winds  the  vertical  velocity  spec¬ 
trum  is  somewhat  steeper  at  high  octaves  and  appears 
to  flatten  out  more  gradually  at  lower  octaves  than  the 
spectra  found  in  weak  mean  wind  conditions. 

Interestingly,  none  of  the  other  spectra  differ  signif- 


9454 


BACMEISTER  ET  AL.:  STRATOSPHERIC  HORIZONTAL  WAVENUMBER  SPECTRA 


a  S  wavenumber(cyc  m“)  )  wavenumberfcyc  m'1) 

10“  10“  10'*  10“  10“  10“ 


0123456789  10  0123456789  10 

Octave  (i/)  Octave  (v) 


Figure  10.  Log-averaged  DWT  power  spectra  of  (a)  ozone  variance  and  (b)  vertical  velocity 
sorted  according  to  magnitude  of  buoyancy  frequency  N.  Two  subsets  of  the  data  are  constructed: 
one  subset  for  measurements  taken  where  N  >  0.023  rad  s-1  according  to  the  NMC  analysis  and 
one  where  N  <  0.020  rad  s-1.  The  high  stratification  group  is  denoted  by  the  solid  triangles, 
and  the  low  stratification  group  is  denoted  by  the  solid  circles. 


icantly  when  sorted  by  mean  horizontal  velocity.  The 
high-mean  velocity  and  low-mean  velocity  power  spec¬ 
tra  for  horizontal  velocity,  potential  temperature  and 
ozone  are  nearly  indistinguishable  from  each  other  or 
from  the  set-averaged  spectra.  Figure  9a  shows  average 
spectra  of  zonal  wind  sorted  according  to  mean  veloc¬ 
ity.  For  horizontal  velocity  and  potential  temperature 
both  the  shapes  and  the  magnitudes  of  the  high-mean 
and  low-mean  velocity  spectra  are  the  same.  For  ozone 
there  is  a  small  increase  in  variance  at  all  horizontal 
scales  for  the  high-mean  velocity  spectra.  However,  the 
shapes  of  the  average  ozone  spectra  for  high-mean  and 
low-mean  velocity  conditions  are  nearly  identical. 


This  result  is  somewhat  at  odds  with  the  findings  of 
Nastrom  et  al.  [1987],  who  obtained  horizontal  veloc¬ 
ity  spectra  with  slopes  steeper  than  -5/3  primarily  over 
rough  terrain  in  high  winds.  In  contrast,  we  see  little 
difference  in  spectral  slope  between  high-  and  low-wind 
environments.  Our  mean  horizontal  velocity  spectra 
under  all  conditions  are  closer  to  the  high-topography, 
high-wind  spectra  of  Nastrom  et  al.  [1987]  than  to  the 
GASP  climatological  spectra.  Although  we  have  not 
isolated  the  oceanic  segments  in  our  data,  we  are  con¬ 
fident  that  they  would  yield  spectra  very  close  to  the 
set-average  spectra  in  Figure  4.  We  make  this  inference 
because  the  slope  histograms  in  Figure  6  are  sharply 
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peaked  despite  the  variety  of  geographical  and  mete¬ 
orological  environments  encountered  during  AASE  II, 
SPADE,  and  ASHOE/MAESA.  In  addition,  over  half 
of  the  total  flight  time  during  the  three  missions  oc¬ 
curred  over  ocean.  Thus  we  conclude  that  typical  con¬ 
ditions  in  the  midstratosphere  are  closer  to  conditions 
of  enhanced  variability  in  the  GASP  data  than  to  cli¬ 
matological  upper  tropospheric  and  lower-stratospheric 
conditions  as  determined  from  the  GASP  data. 

Effect  of  background  Brunt-Vaisala  frequency 

Figure  10  shows  the  effects  of  varying  values  of  the 
inferred  background  Brunt-Vaisala  frequency  N  on  the 
power  spectra  of  ozone  and  vertical  velocity.  The  val¬ 
ues  of  N  used  to  sort  the  spectra  were  derived  from 
NMC  analyses.  This  was  done  to  increase  the  num¬ 
ber  of  days  over  which  N  values  were  available  from 
a  single  source.  The  MTP  was  not  included  on  most 
of  the  SPADE  flights  due  to  payload  weight  consider¬ 
ations.  However,  intercomparisons  of  MTP-measured 
lapse  rates  and  interpolated  NMC  lapse  rates  for  AASE 
II  and  ASHOE/MAESA  show  generally  adequate  agree¬ 
ment  for  1024-s  average  values.  Nevertheless,  this  intro¬ 
duces  a  further  source  of  uncertainty  into  the  analysis. 

The  spectra  were  once  again  divided  into  two  groups. 
One  group  consisted  of  spectra  measured  in  environ¬ 
ments  with  Brunt-Vaisala  frequencies  N  >  0.023  rad  s-1 
and  another  of  spectra  measured  where  N  <  0.020 
rad  s-1.  The  most  striking  changes  associated  with 
stratification  occur  for  the  ozone  power  spectra  in  Fig¬ 
ure  10a.  Mean  variance  for  the  high-stratification  group 
is  reduced  by  nearly  one  decade  at  longer  scales.  At 
shorter  scales  the  difference  in  variance  between  the 
high-  and  low-stratification  groups  is  reduced.  This  im¬ 
plies  a  shallower  ozone  spectrum  in  highly  stratified  en¬ 
vironments.  Power  spectra  of  horizontal  velocity  and 


potential  temperature  (not  shown)  exhibit  only  minor 
sensitivity  to  stratification. 

The  power  spectra  of  vertical  velocity  (Figure  10b) 
show  differences  similar  to  those  seen  in  Figure  9b.  The 
power  spectrum  of  vertical  velocity  in  strongly  strat¬ 
ified  environments  is  similar  to  that  in  environments 
of  low  mean  wind  speed  .  The  likely  explanation  for 
this  lies  in  the  distribution  of  NMC  stratification  fre¬ 
quencies  associated  with  these  missions.  The  histogram 
of  N  for  the  whole  data  set  (Figure  11)  shows  a  dis¬ 
tinct  double-peaked  structure,  with  a  large  peak  at 
N  0.02  rad  s  *  and  a  smaller  secondary  peak  at 
N  ss  0.025  rad  s-1.  The  existence  of  two  distinct  peaks 
reflects  the  relatively  sharp  border  between  tropical  and 
midlatitude  conditions  above  the  subtropical  jet.  NMC 
Brunt-Vaisala  frequencies  higher  than  0.023  rad  s-1  oc¬ 
cur  exclusively  on  tropical  or  equatorward  flights.  The 
relative  sizes  of  the  peaks  simply  reflect  the  higher  num¬ 
ber  of  mid- latitude  and  poleward  flights  during  these 
missions.  Tropical  conditions  tend  to  be  associated  with 
both  weak  horizontal  wind  speeds  and  high  stratifica¬ 
tion.  Thus  there  is  a  large  overlap  between  the  groups 
of  data  intervals  containing  weak  mean  horizontal  wind 
speeds  and  those  containing  strong  mean  stratification. 
The  major  exceptions  to  this  are  the  SPADE  measure¬ 
ments,  and  the  October  1991  to  March  1992  deploy¬ 
ments  during  AASE  II.  These  data  contain  a  high  num¬ 
ber  of  intervals  with  both  light  horizontal  winds  and  low 
Brunt-Vaisala  frequencies. 

To  isolate  the  effect  of  stratification  on  the  power 
spectra  of  vertical  velocity  we  discard  intervals  with 
mean  horizontal  wind  speeds  greater  than  20  m  s"1 
before  grouping  by  stratification.  This  reduces  the  in¬ 
direct  effects  of  mean  horizontal  velocity  on  our  group¬ 
ing  by  stratification.  When  this  is  done  the  differences 


Figure  11.  Histogram  showing  distribution  of  stratification  frequency  values  N  determined  from 
NMC  analyses  for  the  entire  data  set.  The  secondary  peak  of  high  N  values  is  almost  exclusively 
asscociated  with  low-latitude  or  tropical  flight  segments. 
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seen  between  the  high  stratification  and  low  stratifica¬ 
tion  in  Figure  10b  disappear.  Thus  our  velocity  spectra 
do  not  appear  to  vary  significantly  with  mean  strat¬ 
ification,  as  was  inferred  in  earlier  high-wavenumber 
spectral  studies  of  Koprov,  [1965];  Vinnichenko  ei  al., 
[1965];  Myrup,  [1968];  and  Shur,  [1971]  and,  more 
recently,  in  the  analysis  of  GASP  spectra  [Nastrom 
et  al.,  1987].  This  may  be  the  result  of  the  limited 
range  of  Brunt-Vaisala  frequencies  in  our  data.  Al¬ 
most  all  of  our  data  were  taken  in  environments  with 
0.020  rad  s-1  5  N  5  0.025  rad  s-1.  Thus  the  full  range 
of  the  sampled  N  values  may  span  !$20%  of  the  mean. 
The  apparent  sensitivity  of  the  ozone  spectra  to  N  may 
reflect  the  difference  between  tropical  and  midlatitude- 
polar  environments. 

4.  Taylor’s  Hypothesis  in  the  Presence 
of  Gravity  Waves  and  Turbulence 

4.1  Theory 

Traditionally,  aircraft  spectra  such  as  those  used  in 
this  study  have  been  derived  by  implicitly  assuming 
Taylor’s  hypothesis,  so  that  time  variations  measured 
by  the  aircraft  can  be  converted  to  horizontal  length 
scales  by  simply  scaling  the  former  by  the  aircraft’s 
flight  speed  to  obtain  the  latter.  This  requires  that 
the  time  variations  of  the  motions  themselves  are  not 
significant,  or  that  the  aircraft  spans  their  horizontal  ex¬ 
tent  sufficiently  fast  that  they  do  not  vary  temporally. 
Consequently,  aircraft  spectra  are  usually  plotted  as  a 
function  of  horizontal  wavenumber,  as  given  by  a  Taylor 
transformation  of  time  variations  measured  within  the 
aircraft’s  frame  of  reference. 

Here  we  interpret  the  aircraft  data  initially  as  fre¬ 
quency  spectra  measured  from  the  aircraft’s  frame.  This 
allows  us  to  assess  more  accurately  how  gravity  waves 
might  influence  these  spectra,  by  treating  them  as  fre¬ 
quency  spectra  which  have  been  Doppler  shifted  from 
the  intrinsic  (atmosphere-based)  frame  to  the  aircraft 
frame.  The  theory  of  Doppler-shifted  gravity  wave 
spectra  has  been  developed  through  efforts  to  recon¬ 
cile  ground-based  measurements  of  velocity  frequency 
spectra  with  the  estimates  of  both  gravity  wave  and  2- 
D  turbulence  theories,  which  are  typically  formulated 
in  terms  of  intrinsic  frequency  [Schtffltr  and  Liu,  1986; 
Fritts  and  VanZandt,  1987;  Vincent  and  Eckermann, 
1990;  VanZandt  et  al.,  1991;  Gardner  et  al,  1993b]. 
The  2-D  turbulence  theories  always  adopt  Taylor’s  hy¬ 
pothesis  since  this  structure  is  taken  to  be  “frozen”  into 
the  background  flow.  On  the  other  hand,  spectra  pro¬ 
duced  by  gravity  waves  will  generally  include  contribu¬ 
tions  from  motions  which  are  propagating  with  respect 
to  the  background  flow.  If  these  propagation  speeds  are 
comparable  to  the  observing  frame  speed,  then  a  Tay¬ 
lor  transformation  will  not  recover  the  true  horizontal 
wavenumber  spectrum. 

A  general  method  for  calculating  Doppler-shifted 
gravity-wave  spectra  was  given  by  Scheffler  and  Liu 
[1986],  while  Fritts  and  VanZandt  [1987]  presented  more 


illustrative  analytical  solutions  which  incorporated  many 
simplifications  but  retained  the  major  effects.  Gardner 
et  al.  [1993b]  furthered  this  analytical  approach  by  us¬ 
ing  updated  versions  of  the  intrinsic  gravity-wave  spec¬ 
trum  and  by  assuming  isotropic  horizontal  propagation, 
rather  than  the  one-dimensional  (1-D)  horizontal  prop¬ 
agation  considered  by  Fritts  and  VanZandt.  We  adopt 
the  Doppler-shifting  model  of  Gardner  et  al.  [1993b]  in 
our  prediction  of  gravity  wave  frequency  spectra  as  mea¬ 
sured  on  the  aircraft.  Although  their  analytical  model 
assumes  a  hydrostatic  dispersion  relation,  their  results 
differed  little  when  compared  with  more  complete  nu¬ 
merical  calculations  which  used  a  more  accurate  disper¬ 
sion  relation. 

Gardner  et  al.  [1993a, b]  assumed  separability  of  the 
2-D  intrinsic-frequency  and  vertical  wavenumber  spec¬ 
trum  of  horizontal  velocities  produced  by  gravity  waves, 
so  that 

Euv(m,u>)  =  u^otA(m)B(u>), 

where  u2ot  is  the  total  horizontal-velocity  variance. 
Gardner  et  al.  [1993a]  used  a  piecewise-continuous 
model  for  the  vertical  wavenumber  spectrum  A(m), 
shown  with  dashed  lines  in  Figure  12a,  whereas  Gard¬ 
ner  et  al.  [1993b]  used  an  analytical  model  shown  with 
the  solid  curve  in  Figure  12a.  Both  studies  assumed 
B(u>)  <x  w-2,  as  shown  in  Figure  12b.  Similar  models 
to  these  have  been  used  in  recent  spectral  parameteri- 
zations  of  gravity  wave  processes  [Fritts  and  VanZandt, 
1993]  and  agree  broadly  with  spectral  measurements 
of  stratospheric  gravity  waves  [Fritts  and  Chou,  1987; 
Allen  and  Vincent,  1995]  . 

Gardner  et  al.  [1993a]  used  standard  Jacobian  trans¬ 
formations  to  derive  the  1-D  horizontal  wavenumber 
spectrum  Guv(k),  assuming  a  hydrostatic  dispersion  re¬ 
lation  and  isotropic  horizontal  propagation  directions. 
Their  result  is  shown  with  the  dashed  curve  in  Fig¬ 
ure  12c.  Applying  the  hydrostatic  polarization  relations 
to  the  horizontal  velocity  spectrum  yields  the  corre¬ 
sponding  vertical  velocity  spectrum,  i.e. ,  Ew(m,w)  = 
u2N~2Euv(m,u>).  Using  the  transformations  of  Gard¬ 
ner  et  al.  [1993a]  we  derive  the  1-D  horizontal  wavenum¬ 
ber  spectrum  Gw(k)  for  vertical  velocity.  This  is  shown 
with  the  dashed  line  in  Figure  12d.  We  also  derive 
models  for  Guv(k)  and  the  corresponding  vertical  ve¬ 
locity  spectrum  Gw(k)  using  the  models  for  A{m)  and 
B(ui)  adopted  by  Gardner  et  al.  [1993b],  These  are 
plotted  with  solid  curves  in  Figs.  12c  and  12d.  The 
typical  range  of  k  values  resolved  in  these  aircraft  data 
is  shaded. 

Note  that  the  model  horizontal  wavenumber  spectra 
do  not  depend  greatly  on  the  precise  model  used  for 
A{m).  The  Gardner  et  al.  [1993a]  model  for  A(m) 
assumes  ra-5/3  spectra  beyond  mi,  which  is  taken  to 
represent  the  onset  of  inertial  range  3-D  turbulence. 
This  also  yields  k~5^3  spectra  for  k  >  mi,  as  shown  in 
Figure  12c  and  12d,  but  of  course  gravity  wave  theory 
was  not  used  to  derive  this  result.  The  so-called  buoy¬ 
ancy  wavenumber  mj  ‘Ik / Lb  is  discussed  in  more 
depth  below.  It  is  also  worth  noting  that  the  choice  of 
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c):  Guv(k) 


horizontal  wavenumber  k  (eye  m'1) 


b):  B(w) 


intrinsic  frequency  w  (eye  s'1) 

d):  Gw(k) 


horizontal  wavenumber  k  (eye  m'1) 


Figure  12.  Power  spectra  of  horizontal  velocity,  (a)  Power  in  horizontal  velocity  as  a  function 
of  vertical  wavenumber  m.  Power  peaks  near  vertical  wavenumber  m*.  Power  laws  for  spectrum 
away  from  m*  are  indicated  on  the  figure,  b)  Power  in  horizontal  velocity  as  a  function  of  intrinsic 
(atmosphere  based)  frequency  w.  Power  is  assumed  to  be  zero  for  ui  <  /,  where  /  is  the  local 
rotation  frequency,  and  for  w  >  N,  where  N  is  the  local  Brunt- Vaisala  frequency.  Between  f 
and  N,  power  varies  as  u>~2.  (c  and  d)  Doppler-shifted,  Taylor-transformed  power  spectra  for 
(c)  horizontal,  and  (d)  vertical  velocities  as  a  function  of  horizontal  wavenumber.  Shaded  strip 
indicates  scale  range  covered  in  this  study.  The  amplitudes  of  the  model  spectra  of  horizontal 
velocity  are  chosen  arbitrarily  in  this  study. 


the  model  intrinsic-frequency  spectrum  is  directly  re¬ 
flected  in  the  horizontal  wavenumber  spectrum  for  hor¬ 
izontal  and  vertical  velocity  for  /i«  <C  k  m«,  where 
we  have  Guv(k)  oc  k~2 .  Other  choices  for  the  exponent 
in  the  intrinsic  frequency  spectrum  would  change  the 
resulting  horizontal  wavenumber  spectrum  as  well.  In 
other  words  B{u>)  oc  u>~a  will  yield  Guv(k)  oc  k~a  and 
Gw(k)  oc  k~a+ 2  for  <C  &  < C  m*  for  any  choice  of  a. 
The  k~3  behavior  predicted  for  k  >  m*  in  both  horizon¬ 
tal  and  vertical  velocity  power  spectra  results  from  the 
cutoff  of  wave  energy  at  ui  =  N  and  the  m-3  behavior 
assumed  for  energy  at  m  >  m»  [Gardner  et  al.,  1993a]. 

Following  Fritts  and  VanZandt  [1987],  Gardner  et 
al.  [1993b]  found  that  their  Doppler-shifted  frequency 


spectra  were  conveniently  characterized  in  terms  of  a 
Doppler-shifting  parameter 

_  m*(U  -V cos  <j>) 

N  '  W 

where  U  is  the  background  wind  speed,  V  is  the  air¬ 
craft’s  speed,  and  <j>  is  the  azimuth  angle  between  the 
wind  and  aircraft  velocity  vectors.  The  parameter  m* 
is  the  vertical  wavenumber  at  or  near  which  the  hori¬ 
zontal  velocity  power  spectrum  peaks  (see  Figure  12a), 
and  is  generally  predictable  and/or  measurable  (and  rel¬ 
atively  constant)  at  a  given  altitude  [Hines,  1991;  Fritts 
and  VanZandt,  1993;  Allen  and  Vincent,  1995].  Thus 
the  strength  of  the  effects  depends  mainly  on  the  veloc- 


9458 


BACMEISTER  ET  AL.:  STRATOSPHERIC  HORIZONTAL  WAVENUMBER  SPECTRA 


Models  of  Gravity-Wave  Spectra  measured  by  Aircraft 

U  =  mean-wind  velocity:  V=aircraft  velocity 


Doppler**shifted  Frequency  Spectra . 

Q  =  to  +  k  (U  -  V) 
(spectra  of  aircraft  time  series) 

Fuv(Q) 


aircraft  frequency  Q  (eye  s'1) 


k  =  Q/V 


(Taylor  transformed  freq.  spectra) 


Gov(k) 


horizontal  wavenumber  k  (eye  m"1) 


FW(Q) 


aircraft  frequency  Q  (eye  s"1) 


Gw(k) 


horizontal  wavenumber  k  (eye  m'1) 


Figure  13.  Theoretically  derived  horizontal  wavenumber  spectra  for  horizontal  and  vertical 
velocities.  Left  panels  show  frequency  spectra  of  horizontal  velocity  (top  left)  and  vertical  velocity 
(bottom  left)  which  would  be  measured  in  aircraft  traveling  at  different  speeds  with  respect  to 
the  atmosphere.  Thin  solid  line  shows  spectra  for  airspeed  of  0.1  m  s-1  and  should  therefore 
be  close  to  the  unshifted  or  intrinsic  frequency  spectra  shown  in  Figure  15.  As  relative  airspeed 
increases,  Doppler  spreading  smooths  out  the  spectra.  For  airspeeds  above  100  m  s-1  (long 
dashes  and  thick  solid  line),  the  airplane  speed  exceeds  the  propagation  speed  of  a  greater  and 
greater  fraction  of  the  wave  energy  in  the  assumed  spectrum.  Doppler  spreading  ceases  to  affect 
the  shape  of  the  measured  spectra  in  a  significant  way.  Right  panels  show  airplane  measured 
frequency  spectra  transformed  to  horizontal  wavenumber  spectra  using  a  Taylor  transformation. 
Parameters  for  assumed  gravity  spectrum  are  N  —  0.02  rad  s-1  and  m*  =  27r/(3  km).  A  value 
of  /  correpsonding  to  a  latitude  of  40°  is  used. 
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ity  difference  in  (4).  One  should  note  that  Friits  and 
VanZandt  [1987]  and  Gardner  et  al.  [1993b]  considered 
ground-based  frames  exclusively,  and  so  the  frame  speed 
term  V  cos  <j>  does  not  appear  in  their  analysis  since  it 
is  always  zero. 

To  simplify  the  discussion  as  much  as  possible,  we 
shall  consider  the  aircraft  and  wind  velocities  to  be  par¬ 
allel  along  an  east-west  axis,  since  U  is  generally  zonal. 
By  allowing  the  speeds  U  and  V  to  become  positive  or 
negative  zonal  velocities,  then  (4)  can  be  simplified  to 
yield 

*-=4^2.  (5) 

Note  that  when  U  =  V ,  P  =  0  and  no  Doppler  shifting 
occurs,  as  expected. 

The  left  panels  of  Figure  13  show  the  frequency  spec¬ 
tra  of  horizontal  velocities,  Fut)(Cl )  (top),  and  vertical 
velocities,  Fw(£l)  (bottom),  as  measured  on  an  aircraft 
moving  at  seven  different  speeds  with  respect  to  the 
background  wind,  \  U  —  V  \  [Gardner  et  al.,  1993b],  We 
choose  m*  =  27r(3  km)-1  and  N  =  0.02  rad  s-1,  as 
they  are  typical  values  in  the  lower  stratosphere  [Friits 
and  Chou,  1987;  Tsuda  et  al.,  1991;  Allen  and  Vin¬ 
cent,  1995]  see  also  Figure  11.  Note  the  different  fre¬ 
quency  spectral  shapes  that  occur  in  all  seven  cases. 
We  have  computed  these  model  spectra  by  assuming 
equal  amounts  of  wave  energy  propagating  in  all  hor¬ 
izontal  directions  in  the  intrinsic  frame.  This  is  often 
a  good  approximation,  since  the  accumulative  effect  of 
critical  levels  in  filtering  wave  propagation  directions  is 
usually  limited  in  the  lower  stratosphere  (although  not 
always  [see  Eckermann,  1995]).  The  ultimate  Taylor 
transformation  to  a  horizontal  wavenumber  spectrum, 
however,  proves  relatively  insensitive  to  these  details. 

The  corresponding  panels  on  the  right  of  Figure  13 
show  the  result  of  applying  a  Taylor  transformation  to 
the  aircraft  frequency  spectra  on  the  left,  whereupon 
horizontal  wavenumbers  k  are  given  by  k  —  Vl/V ,  and 
the  horizontal  wavenumber  spectrum  is  given  by  [e.g., 
Vincent  and  Eckermann,  1990] 

Guv(k)  =  Fuv(kV)  V,  (6) 

and  likewise  in  transforming  from  Fw(Cl)  to  Gw(k).  If 
we  assume  that  the  variability  in  potential  temperature 
8'  and  trace  gas  mixing  ratios  x  is  due  to  vertical  grav¬ 
ity  wave  displacements  Q'  acting  on  a  mean  vertical  dis¬ 
tribution  of  potential  temperature  0(z)  or  mixing  ratio 
X(z),  i-e., 

then,  using  the  hydrostatic  polarization  relations  we  can 
derive  transformed  variance  spectra  for  potential  tem¬ 
perature,  ozone  and  N2O.  Using  |  C  |=|  (u^^>,)  I  /N, 
these  become 

for  potential  temperature,  and 


GxW=(^x)  Guv(k) 

for  ozone  or  N2O.  Thus  linear  gravity  wave  theory  pre¬ 
dicts  the  same  spectral  shape  for  potential  temperature 
and  ozone  or  N2O  variance  as  it  does  for  horizontal  ve¬ 
locity  variance. 

We  note  from  Figure  13  that  when  |  U  —  V  \  >  30  m  s-1 
(/?~3),  the  frequency  spectra  Taylor-transform  to  over- 
lying  wavenumber  spectra  which  are  identical  to  the 
model  spectra  in  Figs.  12c  and  12d.  This  indicates  that, 
in  these  cases,  the  Taylor  transformation  is  valid  and 
recovers  the  true  gravity  wave  horizontal  wavenumber 
spectrum.  However,  when  |  U  —  V  \  ^  30  m  s-1  (/?<3), 
the  inferred  horizontal  wavenumber  spectra  in  Figure 
13  differ  in  shape.  In  these  cases,  Taylor’s  hypothesis 
is  invalid,  since  the  temporal  oscillations  of  the  waves 
themselves  are  significant  and  contaminate  the  inver¬ 
sion.  These  analytical  results  confirm,  via  a  somewhat 
different  but  essentially  equivalent  route,  the  numerical 
simulations  of  Gardner  and  Gardner  [1993].  Assuming 
U  =  0,  they  concluded  that  aircraft-derived  horizontal 
wavenumber  spectra  produced  by  lower-stratospheric 
gravity  waves  would  not  be  distorted  significantly  when 
|V|»  6  ms-1  (/?  »  1). 

Fortunately  |V|ss  200  ms-1  in  these  aircraft  exper¬ 
iments,  which  usually  ensures  \U  —  V  cos0|jj  50  m  s-1 
when  one  takes  into  account  mean  winds  and  the  cos  <f> 
effects.  Hence  Taylor’s  hypothesis  appears  valid  for 
any  gravity  waves  in  the  data.  Since  the  Taylor  hy¬ 
pothesis  is  always  valid  for  2-D  turbulence,  and  should 
also  hold  for  inertial  and  buoyancy-modified  3-D  tur¬ 
bulence  at  these  large  aircraft  speeds  [e.g,  Powell  and 
Elderkin,  1974;  Wyngaard  and  Clifford,  1977],  then  a 
Taylor  transformation  to  infer  wavenumber  spectra  can 
be  applied  fairly  confidently  to  these  aircraft  data,  re¬ 
gardless  of  their  dynamical  content. 

4.2  Comparison  With  Aircraft  Velocity  Spectra 

The  Taylor-transformed  spectra  on  the  right  of  Fig¬ 
ure  13  faithfully  reproduce  the  model  horizontal 
wavenumber  spectrum  of  gravity  waves  when  \  U  —  V  \ 
is  large.  On  inspecting  these  curves  we  note  that,  at  in¬ 
termediate  horizontal  scales,  Fuv(k )  varies  as  k~2,  while 
l‘\i> (k)  is  essentially  flat.  At  larger  k  both  spectra  ex¬ 
hibit  a  break  towards  k~ 3  shapes.  This  break  occurs 
at  k  =  m,  in  Figure  13  [e.g.,  Gardner  et  al.,  1993a] 
or  at  an  approximate  Fourier  scale  of  2^5  km,  assum¬ 
ing  a  typical  range  of  stratospheric  m*  values.  A  di¬ 
rect  comparison  of  theoretically  derived  aircraft  veloc¬ 
ity  spectra  with  the  observed  aircraft  spectra  is  shown 
in  Figure  14.  The  most  significant  result  is  the  suc¬ 
cessful  prediction  of  the  horizontal  scale  at  which  the 
spectral  index  changes  value,  that  is,  at  a  horizontal 
scale  corresponding  to  k  —  m* .  As  discussed  in  section 
4.1,  this  break  in  the  horizontal  wavenumber  spectra 
is  a  direct  consequence  of  the  cutoff  in  wave  energy  at 
w  =  N,  which  was  assumed  in  the  model  intrinsic  fre¬ 
quency  spectrum.  The  location  of  the  spectral  break 
at  a  given  horizontal  scale  is  determined  by  the  choice 
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of  m,  in  the  model  spectrum.  However,  balloon  and 
profiler  data  consistently  indicate  values  of  m*  between 
2ir(2  km)-1  and  27r(5  km)-1  for  the  lower  stratosphere 
[e.g,  Fritts  and  Chou,  1987;  Tsuda  et  al,  1991;  Allen 
and  Vincent,  1995].  This  range  of  m,  values  leads  to 
breaks  in  the  model  horizontal  wavenumber  spectra, 
which  are  in  good  agreement  with  those  observed  in 
the  ER-2  data. 

The  values  of  the  spectral  indices  within  each  regime 
are  also  well  predicted  by  the  model  spectra.  For  high 
wavenumbers,  Ic  >  m*,  the  observed  spectra  behave  as 
k~2  5 ,  which  is  close  to  the  k~3  behavior  predicted  by 
the  model  gravity  wave  spectra.  The  high  wavenumber 
k~3  behavior  in  the  model  spectra  is  related  to  the  high 
wavenumber  m~3  tail  in  the  assumed  vertical  wavenum¬ 
ber  spectrum  A(m),  which  is  well  established  from  ob¬ 
servations.  For  smaller  horizontal  wavenumbers  the  ob¬ 
served  and  predicted  spectra  are  also  in  good  agree¬ 
ment.  However,  at  these  longer  scales  the  spectral  index 
of  the  model  velocity  spectra  is  directly  related  to  the 
assumed  shape  of  intrinsic  frequency  spectrum  B(u), 
which  must  be  inferred  from  measurements  of  ground- 
based  frequency.  It  may  also  be  worth  noting  that 
the  differences  between  the  theoretical  vertical  velocity 
spectra  for  m,  =  2tr(2.5  km)-1  and  m„  =  2zr(5  km)-1 
in  Figure  14  are  reminiscent  of  the  differences  between 
the  average  vertical  velocity  spectra  for  conditions  of 
strong  horizontal  velocity  and  those  of  weak  horizontal 
velocity  in  Figure  9b. 

In  this  study,  we  regard  the  amplitude  of  the  initial 
model  spectrum  of  horizontal  velocity  Euv(ui,  m )  as  ar¬ 
bitrary.  However,  the  linear  polarization  relations  for 
gravity  waves  still  determine  the  amplitude  of  the  ver¬ 
tical  velocity  spectrum  once  an  amplitude  for  the  hor¬ 
izontal  velocity  spectrum  has  been  chosen.  Thus  the 
gravity  wave  theory  specifies  the  ratio  of  horizontal  ve¬ 
locity  variance  to  variance  of  vertical  velocity  as  a  func¬ 
tion  of  horizontal  wavenumber.  Figure  14c  compares 
theoretical  predictions  of  Guv(k)/Gw(k)  with  measured 
ratios  of  horizontal  and  vertical  velocity  power  spec¬ 
tral  densities.  The  solid  diamonds  show  the  observed 
ratio  of  the  spectral  density  of  total  horizontal  veloc¬ 
ity  to  that  of  vertical  velocity.  As  in  section  3.3,  the 
spectral  density  of  total  horizontal  velocity  is  the  sum 
of  the  zonal  and  meridional  components  (Figure  4a)  at 
each  octave.  The  thin  line  in  Figure  14c  shows  the  ra¬ 
tio  of  the  spectral  density  of  zonal  velocity  variance  to 


that  of  vertical  velocity.  However,  assuming  an  isotropic 
gravity  wave  spectrum,  the  appropriate  comparison  be¬ 
tween  our  theory  and  observations  is  with  the  ratio  of 
total  horizontal  velocity  to  vertical  velocity.  The  the¬ 
ory  predicts  a  ratio  approaching  3.0  for  horizontal  wave 
numbers  k  >  m, ,  and  predicts  rapidly  increasing  ratios 
at  lower  wavenumbers.  In  general,  the  theoretical  ratio 
at  k  >  m,  approaches  (5  —  p)/{ 3  —  p )  for  B(u>)  oc  u>~p . 
The  observed  ratio  of  the  set-averaged  horizontal  veloc¬ 
ity  spectra  agrees  reasonably  well  with  the  prediction  of 
the  gravity  wave  theory  at  longer  scales.  The  agreement 
between  observations  and  theory  seems  to  be  somewhat 
better  for  m,  =  27r(2.5  km)-1,  although  an  even  larger 
value  of  m*  may  be  indicated  by  the  data.  At  smaller 
scales  (<6  km)  the  theoretical  ratio  begins  to  deviate 
significantly  from  the  observed  ratio.  The  theoretical 
ratios  approach  a  value  of  3.0  for  scales  less  than  2n/m *, 
while  the  observed  ratios  are  close  to  6.  Although  dis¬ 
crepancies  occur  at  shorter  scales,  the  general  shape  of 
the  observed  ratio  curve  is  close  to  the  theoretical  pre¬ 
diction.  Overall,  the  agreement  obtained  here  between 
the  theoretical  and  observed  horizontal-to-vertical  ve¬ 
locity  ratio  is  much  closer  than  that  in  previous  studies 
[e.g.,  Gage  et  al.,  1986].  In  their  analysis  of  radar  data, 
Gage  et  al.  [1986]  report  differences  of  up  to  2  orders  of 
magnitude  between  the  observed  ratio  of  horizontal  to 
vertical  velocity  and  that  predicted  by  a  gravity  wave 
theory. 

Several  explanations  are  possible  for  the  discrepancy 
between  observed  and  theoretical  velocity  ratios  in  this 
study.  First,  the  theoretical  ratio  Guv{k)/Gw{k )  in  the 
high-wavenumber  k~ 3  tail  is  somewhat  sensitive  to  the 
assumed  shape  of  B(w).  Steeper  B(u>)  will  lead  to  larger 
Guv(k)/Gw(k)  in  the  tail  portion  of  the  spectrum.  How¬ 
ever,  there  is  no  evidence  for  B(uj)  steeper  than  the  —2 
power  law  assumed  here.  Second,  if  there  are  preferred 
directions  of  propagation,  then  the  total  horizontal  ve¬ 
locity  variance  may  be  closer  to  either  the  individual 
meridional  or  zonal  velocity  variances,  which  would  lead 
to  lower  observed  horizontal-to-vertical  velocity  ratios 
at  high  k.  This  is  evident  in  the  observed  zonal  to  ver¬ 
tical  velocity  ratio  shown  in  Figure  14c.  Finally,  it  is 
possible  and  quite  likely  that  our  gravity  wave  theory 
is  incomplete  at  short  horizontal  scales.  Gravity  waves 
at  these  scales  are  likely  to  be  near  saturation  as  well 
as  being  nonhydrostatic. 

The  possibility  that  quasi-2-D  turbulent  motions  are 


Figure  14.  Comparison  of  measured  log-averaged  DWT  horizontal  wavenumber  spectra  of 
zonal  and  vertical  velocity  (from  Figures  4b  and  6b  respectively)  with  theoretical  predictions  for 
spectra  measured  by  an  airplane  travelling  at  200  ms-1,  (a)  Comparison  for  vertical  velocity. 
Short-long  dashed  curve  shows  prediction  for  m,  =  27t/(2.5  km).  Long  dashes  show  prediction  for 
m«  =  2tt/(5  km),  (b)  Same  as  Figure  14a  except  for  zonal  wind.  It  should  again  be  noted  that  the 
amplitude  of  the  model  spectra  is  arbitrary,  (c)  Ratio  of  horizontal  velocity  variance  to  vertical 
velocity  variance  as  a  function  of  horizontal  wavenumber.  Solid  diamonds  show  measured  ratio 
obtained  from  set-averaged  DWT  power  spectra  of  total  horizontal  velocity  and  vertical  velocity. 
Solid  curve  shows  observed  ratio  of  zonal  velocity  to  vertical  velocity.  Short-long  dashed  curve 
shows  theoretically  predicted  variance  ratio  for  to,  =  27r/(2.5  km)  as  a  function  of  horizontal 
wave  number.  Long  dashed  curve  shows  prediction  for  to,  =  2tt/(5  km). 
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responsible  for  the  apparent  excess  of  horizontal  kinetic 
energy  cannot  be  discounted  [ Gage  et  al,  1986].  How¬ 
ever,  if  this  is  the  case  then  it  must  be  explained  why 
the  excess  appears  most  pronounced  at  scales  shorter 
than  6  km,  and  why  little  excess  of  horizontal  kinetic 
energy  is  apparent  at  scales  longer  than  6  km.  It  is  cu¬ 
rious  that  the  observed  horizontal-to-vertical  velocity 
ratio  and  the  theoretical  ratio  differ  by  a  fairly  uniform 
factor  of  2  at  short  scales.  This  means  that  any  excess 
horizontal  kinetic  energy  must  also  have  a  A~3  spectrum 
at  small  scales,  whatever  its  dynamical  origin.  Further¬ 
more,  this  excess  horizontal  energy  must  also  break  to 
a  shallower  spectrum  at  around  A  =  m* .  Otherwise,  it 
would  quickly  dominate  the  observed  horizontal  veloc¬ 
ity  spectrum,  which  has  a  A-1  5  to  k~ 2  at  scales  longer 
than  k  =  m*. 

Our  measured  velocity  wavenumber  spectra  at  large 
k  are  very  similar  to  those  derived  from  high-resolution 
tropospheric  and  lower  stratospheric  measurements  over 
30  years  ago  by  Shur  [1962],  and  later  by  Vinnichenko 
et  al.  [1965],  Koprov  [1965],  Myrup  [1968],  Shur  [1971], 
Lilly  and  Lester  [1974],  and  Berman  [1976].  In  each  of 
these  studies,  horizontal  wavenumber  spectra  of  veloc¬ 
ity  with  spectral  indices  more  negative  than  -5/3  were 
obtained.  These  observations  were  clearly  inconsistent 
with  the  onset  of  inertial-range  3-D  turbulence.  It  was 
through  efforts  to  explain  these  anomalous  data  that  the 
Lumley-Shur  theory  of  a  buoyancy-modified  subrange 
of  turbulence  first  arose  [SAur,  1962;  Bolgiano ,  1962; 
Lumley,  1964].  This  theory  was  subsequently  refined  in 
several  later  studies  [e.g.,  Weinstock,  1978,  1985b]  and 
has  sometimes  been  invoked  as  a  possible  explanation 
of  m~3  horizontal  velocity  spectra  observed  at  larger 
vertical  wavenumbers  m  [e.g.,  Weinstock,  1985a;  Sidi 
and  Dalaudier ,  1989],  a  spectral  region  more  tradition¬ 
ally  ascribed  to  interacting  and/or  dissipating  gravity 
waves  [e.g.,  Smith  et  al.,  1987;  Weinstock,  1990;  Hines, 
1991;  Gardner,  1994;  Zhu,  1994].  In  fact,  Lumley  [1964] 
argued  that  his  buoyancy-subrange  theory  was  appli¬ 
cable  regardless  of  whether  the  motions  were  gravity 
waves  or  turbulence. 

A  buoyant-subrange  velocity  spectrum  sets  in  at 
wavenumbers  smaller  than  the  buoyancy  wavenumber 
kB  =  2tt/Lb,  given  by 


Lb  = 

cLo’ 

(7) 

To  = 

i j_ y/2 

\N3J  ’ 

(8) 

where  e  is  the  energy  dissipation  rate,  C  is  some  scal¬ 
ing  constant,  and  Lq  is  the  so-called  Ozdimov  scale. 
The  scaling  constant  C  was  not  determined  by  Lum¬ 
ley  [1964].  Weinstock  [1978]  developed  the  buoyant- 
subrange  theory  further,  predicting  kinetic  energy  spec¬ 
tra  ~AP,  where  -4/3  <  p  <  —2.5  at  scales  between  LB 
and  10Tb.  He  also  derived  a  value  C  =  C(um/uo)~ 3/2, 
where  C  «  0.62,  uq  is  the  rms  velocity  within  the  iner¬ 
tial  subrange,  and  um  is  the  rms  velocity  over  the  entire 
spectrum  from  the  scale  of  the  energy  source  to  k  —  oo. 


Weinstock  [1978]  suggested  um/u0S 3  (Tb~50To)  for 
weak  decaying  turbulence.  However,  based  on  sub¬ 
sequent  atmospheric,  oceanic  and  laboratory  studies 
which  indicated  that  these  Lb/Lq  ratios  were  too 
large,  Weinstock  [1992]  derived  a  new  relation 
C  =  (0.16-ftio)-3/4,  where  Rio  is  the  turbulent  gradient 
Richardson  number.  Using  stratospheric  values  for  N , 
and  given  the  possible  wide  variations  in  e,  um/uo  and 
Rio  (and  hence  C),  (7)-(8)  yields  a  range  of  values  for 
the  buoyancy  scale  Lb  ~  10  —  300  m  in  the  strato¬ 
sphere  [e.g.,  Weinstock,  1985a;  Hocking,  1985].  This 
would  be  consistent  with  our  observed  spectral  indices 
of  approximately  —2.5  for  both  horizontal  and  vertical 
velocity  at  scales  below  3.2  km.  Thus  an  ambiguity  in 
interpretation  exists  at  the  high  k  range  of  our  observa¬ 
tions.  The  theory  of  buoyant  subrange  turbulence  and 
gravity  wave  theory  give  similar  results.  However,  this 
agreement  may  reflect  a  true  ambiguity  in  the  nature  of 
the  motion  at  these  scales. 

The  shapes  of  the  aircraft  horizontal  velocity  spec¬ 
tra  at  k  <C  2tt/5  km  resemble  those  found  in  previ¬ 
ous  measurements  in  the  troposphere  and  mesosphere, 
with  spectral  indices  typically  in  the  range  -1.5  to  -2 
[  e.g.,  Nastrom  and  Gage,  1985;  Nastrom  et  al.,  1987; 
Fritts  et  al.,  1989;  Kuion  et  al,  1990;  Hostetler  and 
Gardner,  1994],  There  are  few  published  stratospheric 
measurements  of  this  kind,  however.  Spectra  derived 
from  “HICAT”  data  were  collated  from  various  sources 
and  reviewed  inter  alia  by  Weinstock  [19781.  They  ap¬ 
parently  suggested  spectra  closer  to  A~5'3  at  larger 
scales.  Lilly  and  Lester  [1974]  presented  spectra  of 
along-track  Doppler  wind  velocities  and  potential  tem¬ 
perature  which  suggested  k~3  spectra  at  scales  5 10- 
20  km,  although  their  data  were  inaccurate  at  scales 
<2  km.  Murphy  [1989]  presented  various  ER-2  spec¬ 
tra  but  used  them  more  to  verify  the  accuracy  of  the 
measurements.  Hostetler  and  Gardner  [1994]  presented 
spectra  of  relative  density  fluctuations  but  could  not 
measure  their  spectral  indices  with  any  accuracy.  None 
of  these  earlier  experiments  appear  to  have  noted  the 
break  from  spectra  described  by  A-1  5  or  A-2  to  steeper 
spectra  described  by  A-3.  Lilly  and  Lester  [1974]  also 
presented  a  vertical  velocity  spectrum  inferred  from  the 
potential  temperature  using  gravity  wave  arguments, 
which  forced  the  spectrum  toward  A-1  behavior  at  large 
A.  Thus  it  appears  that  direct  measurements  of  strato¬ 
spheric  vertical  velocity  spectra,  as  presented  here,  have 
not  been  reported  previously. 

Thus  our  observations  tend  to  confirm  the  earliest 
tropospheric  observations  of  horizontal  wavenumber 
spectra,  and  agree  less  with  later  observations  which 
claimed  A-5/3  behavior  within  the  buoyancy  subrange 
and  at  even  larger  scales  [e.g.,  Weinstock,  1978;  Nas¬ 
trom  and  Gage,  1985].  The  spectra  appear  to  be  consis¬ 
tent  with  gravity  waves  and,  at  larger  scales,  also  with 
predictions  of  the  Lumley-Shur-Weinstock  theory  of  a 
buoyant  subrange  of  turbulence  and/or  gravity  waves. 

The  spectral  transformations  involved  in  producing 
Figures  13  and  14  assume  that  the  gravity-wave  disper¬ 
sion  relation  holds  everywhere  [Gardner  et  al.,  1993ab]. 
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The  Doppler-spreading  model  of  Hines  [1991]  argues 
that  the  dispersion  relation  is  invalidly  applied  in  trans¬ 
forming  to  the  portion  of  the  spectrum  which  yields  the 
high  wavenumber  Jfc-3  spectra  in  Figures  12-14,  since 
strong  wave-wave  interactions  make  linear  theory  inap¬ 
plicable  here.  However,  Hines  [1991]  cited  modeling  of 
oceanic  gravity  wave  spectra  by  Allen  and  Joseph  [1989] 
in  which  a  cutoff  horizontal  wavenumber  kc  was  postu¬ 
lated,  beyond  which  dynamical  nonlinearities  yielded 
k~3  spectra  for  both  horizontal  and  vertical  velocities 
measured  in  any  Eulerian  frame.  S.  D.  Eckermann  and 
D.  Broutman  (manuscript  in  preparation,  1996)  applied 
Hines’  Doppler-spreading  ideas  to  show  that  they  pre¬ 
dict  a  cutoff  wavenumber  kc  ft  2m*,  and  hence  at  a 
horizontal  wavelength  ~4-10  km,  similar  to  what  we 
observe  here. 

5,  Discussion 

5.1  Aircraft  Spectra  of  Velocity  and  Potential 
Temperature 

The  good  agreement  between  the  theoretically  pre¬ 
dicted  aircraft  velocity  spectra  and  the  observed  aircraft 
spectra  shown  in  Figure  14  is  a  compelling  argument  for 
a  gravity  wave  interpretation  of  the  variablilty  within 
the  range  of  horizontal  scales  covered  by  our  observa¬ 
tions  (100  km  ^  A  ^  1  km).  The  shapes  of  our  ob¬ 
served  and  theoretical  spectra,  in  particular,  those  for 
vertical  wind,  contrast  sharply  with  the  predictions  of 
a  combined  inverse-cascading,  2-D  and  inertial-range, 
3-D  turbulence  theory,  which  would  have  k~ 5/3  every¬ 
where  for  Guv(k),  and  would  yield  a  similarly-shaped 
spectrum  for  Gw(k)  whenever  the  aircraft  flies  nonpar¬ 
allel  to  surfaces  of  constant  potential  temperature  [e.g., 
Gage,  1979;  Gage  and  Nastrom,  1990].  It  should  also 
be  emphasized  that  a  downscale,  2-D  enstrophy  cascade 
characterized  by  a  kinetic  energy  spectrum  ~k~3  can 
coexist  with  our  model  wave  spectrum  without  funda¬ 
mentally  affecting  our  interpretation  of  the  results.  The 
contribution  of  the  2-D  turbulent  motion  to  kinetic  en¬ 
ergy  would  be  greatest  for  the  horizontal  component 
of  kinetic  energy.  Therefore  its  presence  would  proba¬ 
bly  be  reflected  as  an  excess  of  horizontal  kinetic  energy 
over  that  which  would  be  derived  from  the  vertical  com¬ 
ponent  using  a  pure  wave  interpretation.  However,  the 
ratio  of  horizontal-to-vertical  velocity  variance  in  our 
data  does  not  suggest  the  presence  of  large  amounts  hor¬ 
izontal  kinetic  energy  which  cannot  be  accounted  for  by 
the  wave  interpertration  (Figure  14c).  Thus,  while  an 
enstrophy  cascade  may  be  present  in  the  stratosphere, 
its  contribution  to  kinetic  energy  at  scales  less  than  100 
km  appears  to  be  small  compared  to  that  of  gravity 
wave  motions. 

It  seems  reasonable  to  assume  that  a  transition  to  a 
kinetic  energy  spectrum  dominated  by  2-D  turbulence 
will  occur  at  some  wavenumber.  Our  theoretical  spec¬ 
trum  of  gravity  waves  (Figure  12),  which  is  at  least 
partially  supported  by  observation,  assumes  a  complete 
cutoff  of  gravity  wave  energy  at  an  intrinsic  frequency 


u  =  f.  This  results  in  a  flat  horizontal  wavenum¬ 
ber  spectrum  for  horizontal  velocity  at  scales  longer 
than  27t//i»,  where  h*  =  fm,/N,  for  the  horizontal 
wavenumber  spectra  Guv(k)  in  Figures  13  and  14.  In 
reality,  gravity  wave  energy  may  drop  off  more  grad¬ 
ually  as  ui  approaches  /,  leading  to  a  more  gradual 
decrease  in  the  slope  of  the  horizontal  velocity  spec¬ 
trum.  In  any  case,  the  gravity  wave  interpretation  pre¬ 
dicts  flattening  of  horizontal  velocity  spectrum  at  scales 
~  2 x/h*  ft  1000  km  for  typical  stratospheric  values  of 
N  and  m*  and  midlatitude  values  for  /.  Thus  we  would 
predict  that  somewhere  near  this  scale,  contributions 
to  the  energy  spectrum  from  sources  other  than  grav¬ 
ity  waves  would  begin  to  overwhelm  the  gravity  wave 
spectrum.  Spectra  obtained  from  GASP  data  show  a 
transition  from  a  horizontal  velocity  power  spectrum 
following  a  k~ 5/3  law  to  a  spectrum  ~fc-3  at  a  horizon¬ 
tal  scale  near  1000  km  [Gage  and  Nastrom,  1985].  Gage 
and  Nastrom  1985  argue  that  this  is  the  scale  at  which 
the  inverse  energy  cascade  from  mesoscale  energy  injec¬ 
tions  encounters  the  downscale  enstrophy  cascade  from 
synoptic-scale  energy  injections.  An  alternative  expla¬ 
nation  for  the  change  in  spectral  slope  is  that  the  ob¬ 
served  transition  from  k~ 5^3  behavior  to  k~3  behavior 
at  around  1000  km  reflects  a  transition  from  a  gravity- 
wave-dominated  energy  spectrum  at  smaller  scales  di¬ 
rectly  to  one  dominated  by  enstrophy-cascading  2-D 
turbulence  energized  by  synoptic-scale  inputs.  Simi¬ 
lar  conclusions  were  arrived  at  by  Yuan  and  Hamilton 
[1994]  based  on  nonlinear  numerical  simulations  using 
a  shallow  water  model. 

In  the  2-D  turbulence  view,  the  location  of  the  tran¬ 
sition  from  the  inverse  energy  cascade  to  the  downscale 
enstrophy  cascade  must  be  related  to  the  relative  mag¬ 
nitudes  of  the  synoptic-scale  energy  injections  and  the 
mesocale  energy  injections.  As  the  mesoscale  injection 
decreases  in  magnitude  relative  to  the  synoptic-scale 
injection,  the  transition  from  k~ 5/3  to  k~3  behavior 
would  occur  at  higher  wavenumbers  or  shorter  spatial 
scales.  Specifically,  a  one-decade  decrease  in  the  magni¬ 
tude  of  mesoscale  energy  injections  relative  to  synoptic- 
scale  energy  injections,  should  shift  the  transitional 
wavenumber  up  by  3/4  of  a  decade.  It  is  unlikely  that 
the  relative  magnitude  of  synoptic  and  mesoscale  injec¬ 
tions  is  the  same  at  all  heights  in  the  atmosphere.  In 
particular,  inputs  from  small-scale  and  mesoscale  con¬ 
vection  should  decrease  rapidly  with  altitude  above  the 
tropopause.  Thus  it  would  seem  improbable  that  pure, 
2-D  turbulence  would  produce  transitions  from  k~5^3 
to  k~3  behavior  at  the  same  scale  for  all  altitudes.  On 
the  contrary,  above  the  tropopause,  we  would  expect 
the  transition  to  occur  at  higher  wavenumbers  than  in 
the  troposphere.  The  theoretical  gravity  wave  spectra 
shown  in  Figure  12  predict  more  wave  energy  at  low 
horizontal  wavenumbers  as  m*  and,  consequently,  /i* 
decrease.  The  peak  vertical  wavenumber  m*  is  observed 
to  decrease  with  altitude,  presumably  because  nonlin¬ 
ear  effects  become  important  for  longer  vertical  scales  as 
wave  amplitudes  grow  with  altitude  [Smith  et  ai,  1987; 
Hines,  1991].  Thus  the  combined  wave-enstrophy  cas- 
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cade  picture  would  predict  transitions  from  k~5/3  to 
k~3  behavior  at  lower  horizontal  wavenumbers  (longer 
scales)  for  higher  altitudes. 

A  possible  counterargument  to  the  above  scenario 
would  be  that  a  probable  scale  for  energy  injection  into 
the  middle  atmosphere  exists  at  a  horizontal  scale  of 
around  100  km.  This  source  could  be,  for  example, 
breaking  mountain  waves.  Such  a  source  would  increase 
with  altitude  as  more  and  more  waves  break.  This  en¬ 
ergy  could  go  either  into  a  spectrum  of  secondary  grav¬ 
ity  waves  or  into  the  inverse-cascade  regime  of  2-D  tur¬ 
bulence.  This  could  still  be  consistent  with  a  gravity 
wave  interpretation  of  motion  at  scales  much  smaller 
than  100  km. 

5.2  Aircraft  Spectra  of  Atmospheric  Tracers 

The  ozone  and  N2O  spectra  in  this  study  are  ex¬ 
ceptional  in  that  they  tended  to  conform  to  a  fe-5/3 
power  law.  Similar  results  have  been  inferred  from 
most  previous  spectral  studies  of  tracer-like  chemicals 
at  these  horizontal  scales  [Nastrom  et  al,  1986;  Gage 
and  Nastrom,  1986;  Tjemkes  and  Visser,  1994]  includ¬ 
ing  earlier  analyses  of  ER-2  trace  gas  measurements 
from  AAOE  and  AASE  [Murphy,  1989;  Strahan  and 
Mahlman,  1994].  The  k~5/3  behavior  obtained  for  hor¬ 
izontal  wavenumber  spectra  of  trace  gas  mixing  ratios 
has  typically  been  attributed  to  the  inverse  energy  cas¬ 
cade  of  2-D  turbulence  [Nastrom  et  al,  1986]. 

However,  since  our  spectra  of  horizontal  velocity  and 
potential  temperature  do  not  conform  to  predictions 
for  inverse-cascading  2-D  turbulence,  interpretation  of 
the  contemporaneous  ozone  spectra  in  terms  of  2-D 
turbulence  is  made  difficult.  Studies  employing  high- 
resolution  contour  advection  or  other  air-parcel  advec- 
tion  techniques  [e.g.,  Schoeberl  et  al.,  1993a, b;  Pierce 
et  al,  1994  Plumb  et  al,  1994;  Waugh  et  al,  1994; 
Waugh  and  Plumb,  1994]  have  shown  that  large-scale 
(A  ~  1000  km)  2-D  horizontal  motions  in  the  strato¬ 
sphere  rapidly  produce  highly  distorted  tracer  distri¬ 
butions  with  large  variability  at  much  smaller  scales 
(A  <  100  km)  even  when  the  initial  tracer  field  was 
smooth.  Furthermore,  the  existence  of  these  “vor¬ 
tex  filaments”  can  be  confirmed  through  comparison 
of  these  high-resolution  advective  calculations  with  air¬ 
craft  measurements  of  atmospheric  trace  species  [Schoe¬ 
berl  et  al,  1993b;  Plumb  et  al,  1994;  Waugh  et  al,  1994] 
and  with  satellite  measurements  [Pierce  et  al,  1994].  In 
each  of  these  studies,  the  large-scale  atmospheric  hor¬ 
izontal  velocity  fields  are  obtained  from  global  meteo¬ 
rological  analyses  (e.g.,  NMC/CAC).  These  relatively 
coarse  velocity  fields  are  used  to  advect  large  numbers 
of  purely  passive  trace  particles.  Estimates  of  long-lived 
trace  gas  mixing  ratios  for  each  air  parcel  can  be  esti¬ 
mated,  for  example,  by  using  the  air  parcel’s  initial  po¬ 
tential  vorticity  and  potential  temperature  in  conjunc¬ 
tion  with  large-scale  zonally  averaged  reconstructions  of 
mixing  ratio  [Lait  et  al,  1990]. 

None  of  the  studies  cited  above  derived  spectra  for 
their  calculated  tracer  distributions.  However,  the  pas¬ 
sive  nature  of  the  calculations  makes  upscale  energy 


transfer  impossible  for  scales  smaller  than  the  resolu¬ 
tion  limit  of  the  meteorological  analyses  ^2.5°  x  2.5° 
(inverse  energy  cascading  could  be  implicit  in  the  anal¬ 
yses  at  larger  scales).  Overall,  the  geometry  of  these 
high-resolution  advective  calculations  strongly  resem¬ 
bles  that  in  the  final  stages  of  fully  nonlinear  2-D  cal¬ 
culations  by  McWilliams  [1984].  In  these  calculations 
a  small  number  of  isolated,  coherent,  “hard  core”  vor¬ 
tices  emerge  from  random  noise  initial  conditions.  Be¬ 
tween  these  coherent  vortices  the  flow  is  charcterized 
by  essentially  passive  advection  of  vorticity  filaments. 
Strong  shears  distort  and  stretch  the  filaments,  gener¬ 
ating  smaller  and  smaller  scales  (i.  e.  enstrophy  cas¬ 
cading).  Thus,  it  seems  likely  that  tracer  spectra  de¬ 
rived  from  high-resolution,  advective  calculations  in  the 
stratosphere  would  be  those  associated  with  the  down- 
scale  enstrophy  cascade  of  2-D  turbulence.  Lesieur  and 
Herring  [1985]  have  theoretically  derived  tracer  vari¬ 
ance  spectra  which  vary  as  ife-1  in  this  regime.  The 
shallow  slope  of  the  tracer  variance  spectrum  means 
that  enstrophy- cascading,  quasi-horizontal  turbulence 
can  play  an  important  role  in  producing  tracer  variabil¬ 
ity,  while  at  the  same  time  contributing  little  to  the 
velocity  and  potential  temperature  spectra.  Therefore, 
even  if  the  apparent  excess  of  horizontal  kinetic  energy 
at  small  scales  in  Figure  14c  does  not  represent  quasi¬ 
horizontal  turbulence,  quasi-horizontal  turbulence  may 
still  significantly  perturb  trace  gas  distributions  in  the 
lower  stratosphere. 

The  success  of  parcel  advection  models  in  reproduc¬ 
ing  many  of  the  small-scale  (A  <  100  km)  features  ob¬ 
served  in  aircraft  trace  gas  data  [Plumb  et  al,  1994; 
Waugh  et  al,  1994]  suggests  that  enstrophy-cascading 
2-D  turbulence  is  producing  some  of  the  variability  in¬ 
cluded  in  our  spectra.  However,  our  average  power  spec¬ 
tra  for  ozone  are  clearly  steeper  than  the  k~l  derived 
for  tracers  within  enstrophy-cascading  2-D  turbulence 
by  Lesieur  and  Herring  [1985].  It  is  possible  that  chem¬ 
ical  processes  are  modifying  the  spectrum  of  ozone.  Ex¬ 
pected  lifetimes  for  ozone  in  the  lower  stratosphere  typ¬ 
ically  range  from  several  weeks  to  several  months,  but 
enhanced  destruction  of  ozone  may  occur  if  heteroge¬ 
neous  chemical  processing  is  present.  This  cannot  be 
ruled  out  for  flights  near  the  southern  hemipshere  win¬ 
ter  polar  vortex.  However,  N2O  has  a  sufficiently  long 
photochemical  lifetime  in  the  altitude  range  12-20  km 
that  any  chemical  source  and  sink  effects  on  mesoscale 
variability  can  be  safely  excluded.  Thus  it  appears  un¬ 
likely  that  chemical  processes  play  a  key  role  in  pro¬ 
ducing  the  observed  horizontal  wavenumber  spectrum 
of  ozone  and  N2O. 

We  currently  have  no  satisfactory  explanation  for  the 
observed  spectrum  of  ozone  and  N2O  mixing  ratio  vari¬ 
ance  in  the  EE-2  data.  As  stated  earlier,  the  spectral 
indices  for  ozone  and  N2O  are  consistent  with  the  k~5/3 
spectra  predicted  for  a  tracer  in  the  inverse  energy  cas¬ 
cade  of  2-D  turbulence  [Lesieur  and  Herring,  1985;  Nas¬ 
trom  et  al,  1986].  However,  two  arguments  against  this 
interpretation  have  already  been  implied  and  are  reit¬ 
erated  explicitly  here:  (1)  the  existence  of  the  average 
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horizontal  velocity  and  potential  temperature  spectra 
with  indices  clearly  more  negative  than  -5/3  at  scales 
for  which  a  spectral  index  of  -5/3  is  obtained  for  ozone 
and  N2O;  and  (2)  the  success  achieved  by  other  inves¬ 
tigators  in  reproducing  small-scale  features  in  the  air¬ 
craft  trace  gas  data  using  passively  advected  air  parcels 
in  coarse  resolution  velocity  fields  from  global  analyses. 
Strahan  and  Mahlman  [1994]  offer  several  suggestions 
to  account  for  their  k~ 5'3  N2O  spectra.  They  speculate 
that  an  underlying  k*1  tracer  spectrum  could  be  steep¬ 
ened  by  enhanced  small-scale  diffusion  and  dissipation 
induced  by  cascading  enstrophy,  gravity  wave  breaking, 
or  radiative  damping. 

It  is  clear  that  air  motions  induced  by  gravity  waves 
can  also  directly  induce  variations  in  trace  gas  mixing 
ratios  [e.g.,  Gary,  1989  Danielsen  et  ai,  1991;  Chan  et 
al,  1993;  Teitelbaum  et  ai,  1994;  Reid  et  at.,  1994]. 
The  effect  of  gravity  waves  on  a  given  trace  gas  mixing- 
ratio  will  depend  on  the  local  mixing  ratio  gradients. 
For  short  period  gravity  waves,  vertical  gradients  will 
be  most  important,  but  for  longer-period  inertia  gravity 
waves  both  horizontal  and  vertical  motions  will  be  sig¬ 
nificant  so  that  horizontal  mixing  ratio  gradients  must 
be  considered  [e.g.,  Teitelbaum  et  al.,  1994].  In  the  stud¬ 
ies  cited  above,  gravity  wave  motions  were  assumed  to 
act  on  the  large-scale  climatological  gradients  of  trace 
gas  mixing  ratios  to  produce  variability.  This  should 
lead  to  ozone  and  N2O  spectra  resembling  those  of 
potential  temperature  as  argued  in  section  4.1.  How¬ 
ever,  the  average  potential  temperature  spectrum  de¬ 
rived  here  was  somewhat  steeper  ( k~ 2  to  k~25)  than 
those  observed  for  ozone  and  N2O.  Thus  an  explana¬ 
tion  based  purely  on  gravity  wave  displacements  of  the 
background  tracer  field  is  not  entirely  satisfactory  ei¬ 
ther.  A  more  complete  theory  of  trace  gas  variations 
due  to  gravity  wave  motions  may  have  to  include  the  ef¬ 
fect  of  thin  lamina  or  filaments  of  trace  gases  created  by 
quasi- horizontal,  enstrophy-cascading  2-D  turbulence. 
Schematically,  this  would  appear  as 

Q 

x'  =  C'fo  [x(y<z,t)  + X*(x,y,z,t)]  (9) 

where  is  the  observed  trace  gas  variation,  Q'  is  a  wave- 
induced  vertical  displacement,  x  is  the  climatological 
tracer  distribution  and  \*  represents  tracer  variability 
due  to  isentropic  2-D  turbulent  motions  (e.g.,  the  fil¬ 
aments  simulated  by  high-resolution  parcel  advection 
models).  The  difference  between  (9)  and  earlier  studies 
of  gravity-wave-induced  variability  is  in  the  inclusion  of 
the  x*  term. 

The  potential  effects  of  these  filamentary  or  laminar 
structures  on  gravity-wave-induced  variability  can  be 
appreciated  by  examining  a  typical  polar  ozone  profile 
as  in  Figure  15.  The  profile  in  Figure  15  was  taken  on 
June  30,  1994,  at  Macquarie  Island  (159°E,  55°S).  Mac¬ 
quarie  Island  is  situated  near  the  climatological  edge  of 
the  southern  winter  polar  vortex  and  is  thus  an  ideal 
location  for  examining  the  morphology  of  vortex  ero¬ 
sion.  The  profile  in  Figure  15  shows  numerous  large- 
amplitude  perturbations  in  ozone  mixing  ratios  with 


temperature  (C)  vs  height  (km) 


ozone  partial  pressure  (mPa)  vs  height  (km) 


Figure  15.  (Top)  Temperature  and  (bottom)  ozone 
partial  pressure  as  a  function  of  altitude  over  Macquarie 
Island  on  June  30,  1994. 


very  small  vertical  scales.  Analysis  of  temperature  per¬ 
turbations  for  the  same  balloon  flight  strongly  suggests 
that  these  large  ozone  perturbations  are  not  due  to  grav¬ 
ity  wave  motions. 

A  striking  feature  of  the  profile  in  Figure  15  is  that 
the  vertical  mixing  ratio  gradients  associated  with  the 
thin  lamina  are  actually  much  larger  than  the  mean 
background  gradient.  Therefore  the  total  wave-induced 
tracer  variability  in  (9)  may  be  dominated  by  the  x* 
term.  In  purely  stochastic  terms,  an  important  part 
of  the  total  variablity  would  then  be  the  product  of 
(J  a  process  characterized  by  a  k~2'°  or  k~ 2,5  spectrum 
and  dzX*  ■  Two-dimensional  turbulence  theory  suggests 
that  x*  is  characterized  by  a  k~r  spectrum,  where  r  is 
close  to  1  [Lesieur  and  Herring,  1985].  However,  the  ex¬ 
pected  horizontal  wavenumber  spectrum  of  dz  x*  cannot 
be  inferred  from  strictly  2-D  theory.  Furthermore,  the 
spectral  slope  of  a  product  of  two  stochastic  sequences 
may  depend  on  other  charcteristics  of  each  sequence, 
such  as  intermittency,  and  is  not  easily  predictable. 
Thus  a  prediction  of  gravity-wave-induced  tracer  vari¬ 
ability  appears  to  require  a  better  understanding  of  the 
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3-D  morphology  of  vortex  filaments  or  laminae  than  is 
currently  available. 

The  uncertainty  in  interpreting  stratospheric  tracer 
variability  also  implies  that  caution  should  be  excer- 
cis,ed  in  interpreting  horizontal  wavenumber  spectra  of 
tracers  in  other  regions  of  the  atmosphere.  In  particu¬ 
lar,  horizontal  wavenumber  spectra  derived  from  lidar 
measurements  of  the  mesospheric  sodium  layer  [Kwon 
et  al.,  1990;  Hostetler  and  Gardner,  1994;  Qian  et  al, 
1995;  Gardner  et  al.,  1995]  may  include  structures  anal¬ 
ogous  to  quasi- 2-D  tracer  filaments.  Vertical  profiles 
of  Na  density  in  so-called  sporadic  sodium  layers  [e.g., 
Kwon  et  al.,  1988]  bear  similarities  to  stratospheric 
ozone  profiles  such  as  that  in  Figure  15.  Thus  hor¬ 
izontal  wavenumber  spectra  iderived  from  Na  density 
fluctuations  may  not  be  equal  to  those  of  the  kinetic 
energy  spectra  in  the  mesosphere.  Recent  analyses  by 
Qian  et  al.  [1995]  and  Gardner  et  al.  [1995]  report  ap¬ 
parent  inconsistencies  between  horizontal  wavenumber 
spectra  derived  from  Na  lidar  data  and  those  predicted 
by  gravity  wave  theory  using  the  separable  spectra  of 
Gardner  et  al.  [1993a, b].  In  particular,  they  find  no 
spectral  break  to  steeper  slopes  at  k  =  in  the  hori¬ 
zontal  wavenumber  spectra  of  Na  density  fluctuations. 
Intsead,  the  Na  spectra  of  Qian  et  al.  [1995]  and  Gard¬ 
ner  et  al.  [1995]  look  remarkably  similar  to  our  ozone 
and  N2O  spectra,  exhibiting  slopes  between  -1.3  and 
-1.8  at  horizontal  scales  smaller  than  2ir/mm.  It  is  pos¬ 
sible  that  the  observed  shape  of  these  Na  spectra  reflect 
underlying  quasi-horizontal  motions  rather  than  the  ex¬ 
istence  of  a  nonseparable  gravity  wave  spectrum  as  pro¬ 
posed  by  Gardner  [1994]. 

6.  Summary  and  Conclusions 

We  have  presented  horizontal  wavenumber  spectra 
for  vertical  and  horizontal  velocity,  potential  temper¬ 
ature,  ozone,  and  N2O.  We  used  data  obtained  from  a 
high-altitude  ER-2  aircraft  during  three  aircraft  cam¬ 
paigns  to  examine  stratospheric  chemistry  and  dynam¬ 
ics;  AASE  II,  SPADE,  and  ASHOE/MAESA  (Table  1). 
The  range  of  horizontal  scales  covered  by  our  analy¬ 
sis  extends  approximately  from  400  m  to  102  km.  The 
data  were  taken  at  altitudes  near  20  km  during  73  or  so 
flights.  The  typical  duration  of  each  flight  was  between 
6  and  8  hours  with  a  typical  aircraft  cruise  speed  of 
«400  knots  («200  ms-1).  The  total  length  of  the  ER- 
2  flight  tracks  for  the  three  missions  examined  exceeds 
400,000  km.  Although  the  majority  of  flight  time  ex¬ 
amined  occurred  over  oceans,  a  variety  of  geographical 
and  meteorological  conditions  was  sampled  (Figure  1). 
The  results  of  our  spectral  analysis  using  wavelet  trans¬ 
forms  appeared  to  be  very  robust,  with  both  spectral 
amplitudes  as  well  as  spectral  slopes  between  adjacent 
wavenumbers  exhibiting  sharply  peaked  distributions 
about  their  respective  log-averaged  vaules.  Therefore 
we  are  confident  that  our  results  accurately  represent 
“typical”  conditions  in  the  midstratosphere  (altitudes 
~20  km). 


Our  horizontal  and  vertical  velocity  spectra  as  well 
as  our  potential  temperature  spectra  (section  3.1)  be¬ 
have  approximately  as  fc-2  5  for  scales  smaller  than  6-12 
km.  At  scales  longer  than  12  km  the  horizontal  veloc¬ 
ity  and  potential  temperature  spectra  approach  fc-1  5 
to  k~20  behavior,  while  the  vertical  velocity  spectrum 
becomes  nearly  flat.  This  range  of  behaviors  cannot  cur¬ 
rently  be  explained  by  a  single  turbulence  theory.  The 
k~ 2  5  behavior  at  small  scales  is  approximately  consis¬ 
tent  with  theoretical  predictions  for  buoyant  subrange 
turbulence  [e.g.,  Weinstock,  1978].  At  longer  scales  the 
fc-1  5  to  k~2  0  behavior  is  not  inconsistent  with  inverse- 
cascading  2-D  turbulence  [e.g.  Gage,  1979].  We  propose 
an  alternative  based  on  spectral  transformations  of  an 
observationally  supported  intrinsic  frequency-vertical 
wavenumber  gravity  wave  spectrum  shown  in  Figure  12. 
This  spectrum  is  transformed  to  a  horizontal  wavenum¬ 
ber  spectrum  for  horizontal  and  vertical  velocity  us¬ 
ing  the  linear  gravity  wave  dispersion  relation.  Using 
this  approach,  we  obtain  theoretical  spectra  that  closely 
match  the  observed  spectra  for  horizontal  and  vertical 
wind.  The  spectral  indices  at  long  and  short  scales  are 
reproduced  as  well  as  the  shape  of  the  spectra  as  they 
change  from  the  k~2  5  behavior  at  small  scales  to  k~16 
(horizontal  velocity)  and  k°  (vertical  velocity)  at  the 
longer  scales  in  our  analysis. 

The  spectra  of  ozone  and  N20  which  we  derive  in 
section  3.2  behave  differently  from  those  of  the  meteo¬ 
rological  quantities.  The  spectra  for  these  atmospheric 
tracers  exhibit  behavior  between  k~l  5  and  k~2  0  at 
all  scales.  This  is  not  entirely  consistent  with  a  lin¬ 
ear  gravity  wave  interpretation  which  would  predict  at¬ 
mospheric  tracer  spectra  close  to  those  for  horizontal 
velocity  and  potential  temperature.  We  suggest  that 
this  may  be  due  to  the  presence  of  strong  vertical  gra¬ 
dients  of  mixing  ratio  produced  by  enstrophy-cascading 
2-D  turbulence.  The  spectrum  which  would  result  from 
the  interaction  of  gravity-wave  vertical  displacements 
with  such  structures  in  the  tracer  field  is  as  yet  unpre¬ 
dictable. 

The  analysis  presented  here  is  preliminary.  The  main 
purpose  of  the  present  study  was  to  provide  a  reliable 
characterization  of  typical  mesoscale  variability  in  the 
midstratosphere  and  to  place  it  in  the  context  of  earlier 
observational  and  theoretical  studies.  Although  the  ve¬ 
locity  and  potential  temperature  spectra  derived  here 
strongly  suggest  a  gravity  wave  character  for  much  of 
the  observed  variability  in  the  ER-2  data,  the  spectra 
for  atmospheric  tracers  such  as  ozone  and  N2O  do  not 
appear  to  conform  either  to  a  gravity  wave  interpre¬ 
tation  or  to  a  turbulence  interpretation.  A  complete 
theory  of  mesoscale  variability  in  the  stratosphere  will 
have  to  reconcile  the  wave-like  spectra  of  the  meteoro¬ 
logical  quantities  with  the  anomalous  ozone  and  N2O 
spectra.  Wavelet  transforms  are  a  promising  tool  for 
analysis  of  complex  atmospheric  data  sets  such  as  that 
collected  from  the  ER-2.  We  have  restricted  our  atten¬ 
tion  to  power  spectra  or,  in  other  words,  to  the  sec¬ 
ond  moments  of  wavelet  transforms.  This  yields  in¬ 
formation  equivalent  to  that  obtained  from  traditional 
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Fourier  techniques.  However,  higher  order  moments  of 
wavelet  transforms  can  be  used  to  quantify  other  char¬ 
acteristics  of  data  such  as  intermittency  [e.g.,  Katul  et 
al.  1994].  Such  analyses  of  the  ER-2  data  may  shed 
more  light  on  the  dynamical  processes  producing  strato¬ 
spheric  mesoscale  variability. 
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